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We present a measurement of the mass dependence of the forward-backward charge asymmetry 
(Aps) for e'*'e~ pairs produced via an intermediate Z/^* with mass Mee > 40 GeV/c^. We study 
the constraints on the Z-quark coupUngs imposed by our measurement. We analyze an integrated 
luminosity of 72 pb~^ collected by the CDF II detector in pp collisions at -^s — 1.96 TeV at 
the Fermilab Tevatron. A comparison of the uncorrected Afb between data and Standard Model 
Monte Carlo gives good agreement with a x^/DOF of 15.7/15. The couplings measurements are 
also consistent with Standard Model predictions. 

PACS numbers: 13.85.Qk 12.38.Qk 12.15.Ji 12.15.Mm 



I. INTRODUCTION 

The reaction pp £'^£~ , where £ is an isolated high- 
Pj, lepton, is mediated primarily by virtual photons at 
low values of dilepton invariant mass (M^+^-) pri- 
marily by the Z at Mg+g- ~ Mz, and by a combina- 
tion of photons and Z bosons outside these regions. The 
presence of both vector and axial- vector couplings of elec- 
troweak bosons to fermions in the process qq — > Z/j* — > 
£~^£~ gives rise to an asymmetry in the polar angle of the 
lepton momentum relative to the incoming quark mo- 
mentum in the rest frame of the lepton pair. 

At tree level, the process qq — > £'^£~ proceeds via 



an s-channel exchange of either a virtual photon or a 
Z boson. The neutral current coupling of a fermion / 
to the Z boson has vector and axial- vector components: 

J^-'' ~ f{9v + 9Al5)f, where and g{ are the vector 
and axial- vector couplings of the fermion to the Z respec- 
tively. The coupling of the same fermion to the photon is 
purely a vector coupling and its strength is proportional 
to the charge of the fermion Qf. The differential cross 
section for qq £^£^ is obtained by squaring the matrix 
element, integrating over the azimuthal angle, averaging 
over the polarization of the incoming particles, and sum- 
ming over the spin and polarization of the final state 
particles: 
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da{qq -> £+£- 
dcos 9 



+ 



4.9l5i cos ( 
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where C is the color factor, 9 is the emission angle of the 
lepton (anti-lepton) relative to the quark (anti-quark) in 
the rest frame of the lepton pair, and 



1 



cos^ 9w sin 9w s 
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The first and the third terms in Eq. correspond to 
the pure 7* and Z exchange respectively while the sec- 
ond term corresponds to the Z /^* interference. The an- 
gular dependence of the various terms is either cos6' or 
( 1 -I- cos^ 9) . The cos 9 terms integrate to zero in the total 
cross section but induce the forward-backward asymme- 
try. 
Let 
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R 



AA 



QiQ; + 2Q,Q,glg'yRe{x{s)) 
+ 9'v{9l'+9f)\x{s)? 
^ 9'A{9V+9f)\x{s)\^ 
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i?^^ ^ lg\g'A{Q^Q,Re{x{s)) + 2glg'y\x{s)\').{h) 

Then the differential cross section is reduced to the fol- 
lowing simple expression: 

da ^ c'^'^ ^ A {I ^ cos" e)+ApB cos 9] (6) 
a cos 9 3 s 8 

where Rf = R^^ + R^^ and ApB = The meaning 

of the quantities Rf and ApB can be clearly seen. In- 
tegrating Eq. © over cos0, the first term in the square 
brackets integrates to unity, the second integrates to 0. 
Therefore atotai = CRfa^^^ , where <Tq^^ is the total 
QED cross section (the cross section if the Z^ exchange 
amplitude were absent). The quantity Aps can be writ- 
ten as 
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FB 



-1 da 



idcosO-\- L ^ -T^^d cos 9 

JO a cos tf 



dcosi 



TT^d cos 9 

J — 1 a cosy 



ap — ub 

ap + <TB 



(7) 



and is identified as the forward-backward asymmetry, 
where is the number of forward (cos^ > 0) events 
and is the number of backward (cos 9 < 0) events. 

A measurement of ApB can constrain the properties 
of any additional non-Standard Model amplitudes con- 
tributing to qq £+£- [3, and is complementary to 



direct searches for non-Standard Model amplitudes that 
look for an excess in the total cross section. This is par- 
ticularly interesting for Mee above LEP II energies, where 
the measurement is unique to the Tevatron. 

In addition, Eq. ^ shows that depending on the in- 
variant mass, a different combination of vector and axial- 
vector couplings contribute to the differential cross sec- 
tion. Consequently, Ap ^ is a direct probe of the relative 
strengths of the coupling constants between the Z boson 
and the quarks. The invariant-mass dependence of ApB 
is also sensitive to u and d quarks separately, unlike other 
precise measurements of light quark Z couplings in i^N 
scattering f3| and atomic parity violation 5] on heavy 
nuclei. 

In this paper, a number of different comparisons be- 
tween the data and Standard Model expectations are 
presented. The uncorrected M^e, cos 9*, and Aps data 
distributions are compared with the output of our Monte 
Carlo and detector simulation. The first principal result 
is a measurement of Aps in 15 M^e bins using an un- 
folding analysis that doesn't assume a prior Standard 
Model ApB distribution. The second principal result 
is a measurement of three sets of parameters: the Z- 
quark couplings, the Z-electron couplings, and sm^9w- 
In making each of these three measurements, the other 
parameters are held fixed with the values given by the 
Standard Model. The measured Z-quark couplings are 
then used to determine experimental correction factors 
for acceptance and efficiency of dielectron events. The 
correction factors are used for the measurement of ApB 
with Standard Model assumptions. 

The previous measurement from the Collider Detector 
at Fermilab (CDF) 1^ was made with the data taken be- 
tween 1992 and 1995 using the Run I detector. The Run 
I measurement assumed a Standard Model Ap b distribu- 
tion for calculating efficiency and acceptance experimen- 
tal correction factors. The present measurement uses the 
dielectron data taken between March 2002 and January 
2003 with the CDF II detector, corresponding to an in- 
tegrated luminosity of 72 pb^^. 

The paper is structured as follows. A description of 
the detector and an overview of the analysis are given 
in Sec. |n] Event selection and candidate events are dis- 
cussed in Sec. IIIII The estimation and characteristics of 
the backgrounds are described in Sec. IIVI The accep- 
tance and corrections for detector effects are described 
in Sec. The systematic uncertainties are summarized 
in Sec. I VII Finally, the results of the forward-backward 
asymmetry and coupling measurements are presented in 
Sec.rVIll 
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II. OVERVIEW 

This section begins with a discussion of aspects of the 
detector, triggers and data samples that are relevant to 
this measurement. The nature of Z j^* e~^e~ events 
in a hadron collider and the overall strategy of the anal- 
ysis are presented. This paper uses a cylindrical coordi- 
nate system, with the positive z axis oriented along the 
beamline in the direction of the proton's momentum. 



A. Detector and triggers 

The Collider Detector at Fermilab (CDF II) is a 
general-purpose detector designed to study the physics 
of collisions at ^/s = 1.96 TeV at the Fermilab Teva- 
tron collider. Like most detectors used at high-energy 
colliders, it has a cylindrical geometry with axial and 
forward-backward symmetry. A diagram of the inner 
part of the CDF II detector is shown in Fig.^ The inner- 
most part of the detector contains an integrated track- 
ing system with a silicon detector and an open-cell drift 
chamber. A solenoidal magnet surrounding the tracking 
chambers provides a 1.4 T field aligned with the proton 
beam axis. The integrated tracking system is surrounded 
by calorimeters which cover 2tt in azimuth and from —3.6 
to 3.6 in pseudorapidity, rjdet (see Fig.^. Outside of the 
calorimeters is a muon system with coverage from —1.5 
to 1.5 in rjdet- The CDF II detector is a major upgrade to 
the detector that took data until 1996. The entire track- 
ing system subtending \ridet \ < 2 and the plug calorimeter 
subtending 1.1 < jrydetl < 3.6 have been replaced to han- 
dle the higher rate of collisions and increase the capabili- 
ties for physics analyses in Run II. This analysis uses the 
open-cell drift chamber called the central outer tracker 
(COT) and the calorimeters. A more detailed detector 
description can be found in Refs. 0,13, and a description 
of the upgraded detector can be found in Ref. 

The COT detector ^3 is a 96 layer, 3.2 m long open- 
cell drift chamber which, combined with the solenoid, is 
used to measure the momenta of charged particles with 
\Vdet\ < 1. The detector extends from a radius of 40 cm 
to a radius of 137 cm. The 96 layers are divided into 
8 "super-layers", which alternate between super-layers 
where the wires are axial (i.e., parallel to the z axis) and 
super-layers where the wires have a ±2° stereo angle, 
providing three-dimensional tracking. The calorimeter 
consists of a lead-scintillator electromagnetic (EM) com- 
partment with shower position detection backed by an 
iron-scintillator hadronic compartment. The calorime- 
ters are segmented in projective r]det~4' towers pointing 
to the nominal interaction point, at z = 0. While the 
central calorimeter (|»7d£t| < is retained mostly un- 
changed from Run I Ul El El El El El: the plug 
calorimeter with 1.1 < \r]ciet \ < 3.6 is a major compo- 
nent of the Run II upgrade (Fig.[2Jl, and largely follows 
the design of the central detector. Since the calorime- 
ter segmentation is coarse compared to the dimensions 



7,= 1 



HAD CENTRAL CALOIUMETEn 



I CENTRAL CALORIMETER 



IsqlehqidI 




SILCON LAYERS 



FIG. 1: One quadrant of the CDF II tracking and calorimetric 
detectors. The detectors have axial and reflective symmetry 
about z—0. CDF uses a cylindrical coordinate system with the 
z (longitudinal) axis along the proton-beam direction; r is the 
transverse coordinate, and (j) is the azimuthal angle. The de- 
tector pseudorapidity is defined as rjdet = — ln(tan^^), where 
6det is the polar angle relative to the proton-beam direction 
measured from z = Q. The event pseudorapidity is defined as 
rjevt = — ln(tan^^), where O^vt is the polar angle measured 
from the nominal pp collision point in z. The transverse mo- 
mentum (pj,) and energy (Et) are the components projected 
onto the plane perpendicular to the beam axis (p^ = p ■ sin 8 ; 
Et = E-sin6). The missing transverse energy, I^t, is defined 
as the magnitude of —T,iEThi, where hi is a unit vector in 
the perpendicular plane that points from the beamline to the 
ith calorimeter tower. 



of an electron shower, position detectors (shower maxi- 
mum detectors, CES in the central region and PES in 
the plug region) are placed at a depth of approximately 
six radiation lengths, roughly the position of the shower 
maximum, inside the EM calorimeters. These detectors 
measure the position and profile of the showers and help 
differentiate electrons from hadrons. 

The trigger system has undergone a complete redesign 
as a result of the accelerator and detector upgrades. The 
CDF trigger is a three-level system that selects events 
out of a 2.5 MHz crossing rate to be written to magnetic 
tape at a rate of ~75 Hz. The first two levels [3 are 
composed of custom electronics with data paths separate 
from the data acquisition system. The third level 
receives the complete detector information from the data 
acquisition system and runs a version of the reconstruc- 
tion software optimized for speed on a farm of commer- 
cial computers. The Z/j* events used in this measure- 
ment are selected by high-i?T central electron triggers. 
At Level 1, electrons are selected by the presence of an 
energy deposition of Et > 8 GeV in an EM calorime- 
ter tower and a matching two-dimensional (r-(f) plane) 
track with pT > 8 GeV/c reconstructed in the COT by 
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• W ^ ev sample: A sample of 38,000 W ^ ev 
candidates, where the electron is reconstructed in 
the central calorimeter, is used to study the mate- 
rial in the tracking volume, to calibrate the relative 
calorimeter response within a central tower, and to 
check charge biases in measuring electrons. 

• Inclusive electron sample: A sample of 3 million 
central electron candidates with Et > 8 GeV is 
used to calibrate the relative response of the central 
EM calorimeter towers. 




AS = 2.7" 

An = ae / sin e 



FIG. 2: Cross section of upper part of plug calorimeter (top), 
and transverse segmentation, stiowing physical and trigger 
towers in a 30° cf) section (bottom). The logical segmenta- 
tion for clustering purposes is the same except in the outer 
two rings (9 > 30°^, where two neighboring (in azimuth) 7.5° 
towers are merged to match the 15° segmentation of the cen- 
tral and wall calorimeters behind them. 



• Dijet samples: A sample of 1 million dijet events 
(events with at least two jets, each with Et > 20 
GeV) is used to measure the rate at which jets 
fake an electron signature and to estimate the di- 
jet background. A jet is defined as a cluster of 
energy reconstructed in the calorimeter. Triggers 
for the sample require a calorimeter tower with 
Et > 5 GeV at Level 1, a calorimeter cluster with 
Et > 15 GeV at Level 2, and a reconstructed jet 
with Et > 20 GeV at Level 3. Due to the high 
cross section, only 1 in roughly 400 events on aver- 
age were randomly selected to be recorded for this 
trigger. Jet samples with higher trigger thresholds 
(50 GeV, 70 GeV, and 100 GeV at Level 3) are also 
used to cross-check the fake rate for a trigger bias. 



C. Monte Carlo samples 



the extremely Fast Tracker (XFT) "201. At Level 2, the 
electron energy is reconstructed in a cluster of EM tow- 
ers and is required to have Et > 16 GeV. At Level 3, 
a reconstructed EM cluster with Et > 18 GeV and a 
matching three-dimensional track with pT > 9 GeV/c 
are required. At each level, the energy in the hadronic 
towers just behind the EM tower or cluster is required 
to be less than 12.5% of the EM energy. Another trigger 
path that does not require a hadronic energy fraction is 
added to this measurement in order to improve the trig- 
ger efficiency for central electrons with very high Et. The 
integrated luminosity of the data sample for this analysis 
is (72 ±4) pb-i. 

B. Data samples 

Four data samples are employed in this analysis. These 
are described briefly below and in more detail in subse- 
quent sections. 

• Z/'y* — > e~^e~' sample: A sample of 5,200 dielectron 
candidates is used to measure ApB^ calibrate the 
energy scale and resolution of the EM calorimeter, 
and study the material in the tracking volume. 



Monte Carlo generation and detector simulation are 
used to measure the acceptance for the Drell-Yan process, 
model the effect of radiation and detector resolution, de- 
termine the characteristics and amount of background in 
the data sample, and understand systematic uncertain- 
ties on the ApB measurement. PYTHIA |21| and HER- 
WIG ^ generators with CTEQ5L PDFs |23| are used 
for most of the samples. These generate processes at lead- 
ing order and incorporate initial- and final-state QCD 
and QED radiation via their parton shower algorithms 
(HERWIG does not include final-state QED radiation). 
PYTHIA is tuned so that the underlying event (the re- 
maining pp fragments from the collision) and the spec- 
trum of Z bosons agree with the CDF data Two 
matrix element generators, WGAMMA and ALP- 
GEN ,261, are used to check the W -\- X ^ ev -\- X back- 
ground estimate from PYTHIA, where X is a photon or 
hadronic jet. The generator WGAMMA calculates the 
cross section of the pp — > W-f process. It uses electroweak 
helicity amplitudes for production and radiative W 
boson decays, including all interference terms. ALPGEN 
performs the calculation of the matrix element for the 
production of PI^-Kquark and H^-|-gluon final states. The 
detector simulation models the decay of generated par- 
ticles and their interactions with the various elements of 
the CDF detector. A fuU GEANT3 simulation is 
used to simulate the tracking volume. A parameterized 



simulation is used for the calorimeters and muon cham- 
bers. Comparisons between the data and Monte Carlo 
simulation are discussed in Sec. IIII Dl 

There are nine Monte Carlo samples used in this anal- 
ysis, which are briefly described below. 

• Z 1^* — > e+e^ sample: A sample generated with 
PYTHIA is used to calculate corrections due to ac- 
ceptance, bremsstrahlung, and energy resolution, 
and to estimate the systematic uncertainties due to 
the energy scale and resolution. A quarter of the 
sample was generated with M^f-y^ > 105 GeV/c^ to 
reduce the statistical uncertainties associated with 
the Monte Carlo sample in the high-mass region. 

• Z/"f* e^e~ sample for material systematics: 
Three PYTHIA samples are used to estimate the 
change in the measured ApB between the default 
simulation and a modified simulation which adds or 
subtracts 1.5% radiation length (Xq) of copper in 
a cylinder in the central region and ^Xq of iron on 
the face of the plug calorimeter. QCD fragmenta- 
tion is turned off for these samples in order to save 
CPU time. 

• Z/-f* — > T+T^ sample: A PYTHIA sample is used 
to estimate the background due to Z/j* t^t~ . 
TAUOLA m is used to decay r's. 

• Dijet sample: A PYTHIA sample with aU 2 — > 2 
processes is used to understand the characteristics 
of the dijet background. A lower limit of Pj. > 
18 GeV on the transverse momentum in the rest 
frame of the hard interaction is applied. 

• ti sample: A HERWIG sample is used to estimate 
the background due to ti production. 

• Diboson samples: A sample with WW production 
and a sample with WZ production are generated 
using PYTHIA and used to estimate the diboson 
backgrounds. 

• W ^ ev sample: A PYTHIA sample is used to 
estimate the background due to the inclusive W 
production. 

• W^ ei/j sample: A WGAMMA sample is used 
to cross-check the background due to -I- 7 pro- 
duction. 

9 W + q/g ev + q/g sample: An ALPGEN sample 
is used to estimate the background due to -I- 
quark or gluon production. 



D. Strategy of analysis 

This analysis focuses on e'^e' pair production via an 
intermediate Z/'y*. The goal of this analysis is to mea- 
sure ApB as a function of the invariant mass of the di- 
electron pair. The dielectron sample is chosen because of 




FIG. 3: Diagram of pp Z/j*, where one of the quarks 
radiates a gluon or photon imparting transverse momentum 
to the quark. Once the quarks annihilate to Z/'y* , the source 
of the transverse momentum is ambiguous. 



the low backgrounds and the good polar angle coverage 
of electrons in the CDF II detector. 

When the incoming quarks participating in the Drell- 
Yan process have no transverse momentum relative to 
their parent baryons, 9 in Eq. © is determined unam- 
biguously from the four-momenta of the electrons by cal- 
culating the angle that the electron makes with the pro- 
ton beam in the center-of-mass frame of the electron- 
positron pair. When either of the incoming quarks has 
significant transverse momentum, however, there exists 
an ambiguity in the four- momenta of the incoming quarks 
in the frame of the electron-positron pair, since one can- 
not determine the four-momenta of the quark and anti- 
quark individually (see Fig. O). The CoUins-Sopcr for- 
malism [23 is adopted to minimize the effects of the 
transverse momentum of the incoming quarks in Eq. 10 . 
(The magnitude of the effect in the CoUins-Soper frame 
is discussed in Sec. IVII Bl ) In this formalism, the polar 
axis is defined as the bisector of the proton beam momen- 
tum and the negative of the anti-proton beam momentum 
when they are boosted into the center-of-mass frame of 
the electron-positron pair. The variable d* is defined as 
the angle between the electron and the polar axis. Let 
Q (Qt) be the four-momentum (transverse momentum) 
of the electron-positron pair. Pi be the four-momentum 
of the electron, and P2 be the four-momentum of the 
positron, all measured in the lab frame. Then cos6'* is 
given by 



COS! 



VOW' 



iP+P,-p-p+), (8) 



where P^ — -^{Pi ^Pi)i and and P^ represent the 
energy and the longitudinal component of the momen- 
tum, respectively. Forward events are defined as having 
cos 6** > and backward events are defined as having 
cose* < 0. 

The ApB measurement is made in 15 bins of M^e be- 
tween 40 GeV/c^ and 600 GeV/c^. The bin sizes are cho- 
sen based on the detector Mee resolution and the relative 
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Drell-Yan cross sections in each bin. Once the event se- 
lection has been made, the estimated number of forward 
and backward background events are subtracted from the 
candidate sample in each Mf.e bin. The raw forward- 
backward asymmetry (A™^ ) calculated by applying 
Eq. lO to the background subtracted sample. 

The goal is to connect what we measure in the detec- 
tor after background subtraction (A™^) to theoretical 
calculations {A^p^^). If the detector simulation is avail- 
able, a direct comparison can be made between the uncor- 
rected data and a model prediction. These comparisons 
are made in Sec. IVII Al If the simulation is not avail- 
able, the true A^pj^' must be unfolded from A',?^ . A"^^ 
must be corrected for detector acceptances and efficien- 
cies (which sculpt the cos 6* distribution) and for smear- 
ing effects (which move events between M^e bins) to ob- 
tain A^pg^ , which can be compared directly to theoretical 
calculations. The Mee bin migration near the Z pole is 
not negligible, which makes it difficult to unfold A^g"^ in 
those bins. Performing an unconstrained unfolding anal- 
ysis in those bins results in large uncertainties, while ap- 
plying Standard Model or other constraints can bias the 
result. Two different unfolding analyses are performed 
to reconstruct A^^^ with different constraints. Both un- 



folding analyses assume the Standard Model 



da 



to be 



able to calculate the event migration effects. The uncer- 
tainties associated with systematic effects are estimated 
by varying the magnitude of these effects in the Monte 
Carlo simulation and recalculating the results for each 
analysis. 

The first and least constrained analysis is performed by 
fitting for A^g'' with a smoothing constraint. The accep- 
tance and event migration are parameterized (Sec. IV D(l 
to transform Af^' into A-^ 



1 raw 



{AWS)^ = .9(A,z), (9) 
where A = {(<J^)i, (A^/j^)^, . . . , (<r)i5}-(10) 

The probability, for a given Aps, to find a number of 
forward and backward events is the binomial probability 
{P[Afb,N^,N^]). The maximum likelihood method is 
used to obtain the 15 values of A^pg'^ or A fSec. IVIlT^ . 
where the negative log likelihood with a smoothing con- 
straint for bins near the Z pole is defined as 

15 9 

a = Y. -lo9{P[9{A, z), Nf,Nf]) A • ^ S[A.l (11) 



perform a fit between the smeared theoretical calcu- 
lations and A'^pg to extract the couplings. This analysis 
is described in detail in Sec. IVII Dl 

Continuing with the second analysis, Aps can be mea- 
sured by assuming the Standard Model Aps shape with 
the measured Z-quark couplings. In this method, ac- 
ceptance correction factors (ajg,.) are used to correct the 
number of forward and backward events in each M^e bin. 
The correction factor accounts for the effects of detector 
acceptance and efficiency, and bin migration (Sec. IV^ . 
These correction factors depend on the input A^g'' that 
is used to calculate them. To constrain the Standard 
Model bias, the input A^g'^ is constrained to the predic- 
tion with uncertainties from the Z-quark coupling analy- 
sis. Then A^pg'^ can be calculated directly using the cor- 
rection factors fSec. IVIIE)) . Equation {7)) can be rewrit- 
ten: 



A 



tda_\ + 
\dM )i 



FB 



lda_Y 
\ dM ) i 



lda_y 
\dM li 



fda_Y 
\dM li 



where 



\dM ) ^ AMiL{a: 



± 1 

cor ) 



(13) 



(14) 



In Eq. p4l) . AMi is the size of the i-th mass bin, and L is 
the integrated luminosity. N~^{N~) and B^{B~) are the 
number of candidate and estimated background events in 
the forward (backward) region, respectively. Canceling 
common factors, the forward-backward asymmetry can 
be written as 



A 



FB 



TV" 



{aaaT)i 



N~ 



-B~ 



(15) 



1^ 



(acor)i 



The analysis using correction factors has smaller uncer- 
tainties than the unconstrained method, but it assumes 
no non-standard physics outside of deviations in Z-quark 
couplings. This analysis technique is similar to the Run 
I measurement. 

While these analyses make different constraints on the 
form of ApB, it is important to note that they should 
give similar results in bins where there is a negligible 
amount of event migration between M^e bins. At high 
mass, where the bin sizes are very large compared to the 
detector resolution, the results are independent of the 
unfolding method. 



where A, = {(<J^).-i, (<J0„ (<r),+i}, (12) 

S'[Aj] is a regularization function, and A is the regular- 
ization parameter. This analysis makes no assumptions 
about the shape of ApB (aside from it being smooth near 
the Z pole), but it has the largest uncertainties. 

In a second analysis, the parameterized acceptance 
and event migration are also used for the measurement 
of the Z-quark and Z-electron coupling constants, and 
sin^ . We vary the couplings at generator level and 



III. ELECTRON MEASUREMENT 

This analysis requires two electrons (e+,e^) in the 
event, one in the central region, and the other in the 
central or plug region. This section describes the identi- 
fication of central and plug electrons, the event selection 
criteria, the electron energy scale and resolution, and the 
charge identification of electrons. In most cases, discus- 
sion of electrons refer to both electrons and positrons. 
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A. Central electron identification 



Electron identification in the central calorimeter is al- 
most identical to the algorithm used in Run I, since the 
calorimeter is unchanged and the new drift chamber has 
a geometry very similar to the previous one. For a more 
detailed description of the central electron reconstruction 
variables see Ref. ^. An electron candidate is recon- 
structed if there is a central tower with Et > 2 GcV 
and a charged drift chamber track that extrapolates to 
the tower. The adjacent towers on either side in rjdet 
are added to the cluster, and the cluster is not accepted 
if the energy in the hadronic part is more than 12.5% 
of the energy in the EM part. An electron is consid- 
ered within the fiducial region of the detector if its track 
points within 60 cm in z of the center of the detector 
and extrapolates to the calorimeter away from any wedge 
boundaries. (A wedge consists of those towers with the 
same value of 0.) The polar range of electrons in the 
central region is \ridet\ < 1-0. The energy of the electron 
is determined by the total energy it deposits in the EM 
calorimeter. The momentum (pj.) of the electron is de- 
termined by the highest-p^, COT track associated with 
the EM cluster. The track is constrained to the position 
of the beamline in The direction of the electron's 

momentum is taken from the direction of the track and 
is used in the calculations of the transverse component 
of the energy (Et) and the invariant mass of the elec- 
tron pairs. The charge of the electron (q) is determined 
from the curvature of the track. The variables that are 
used to discriminate electrons from hadrons are: (1) the 
ratio of the hadronic energy to the electromagnetic en- 
ergy {Ehad/ Eem)'i (2) the total transverse energy within a 
radius of 0.4 in Ai? = ^Vevt + ^'/'^ of the cluster cen- 
troid, excluding the cluster energy itself (i?^°);(3) the 
ratio of the calorimeter energy to the momentum of the 
track (E/p); (4) the comparison of the lateral sharing of 
energy among the calorimeter towers with that of test 
beam electron data [Lshr]] (6) a comparison of the 
shower profile measured by the shower maximum detec- 
tor with the shower profile measured from test beam elec- 
trons (Xstrip); (5) the distance in r-4> and z between the 
electron shower position measured by the shower maxi- 
mum detector and the extrapolated track position (q • Ax 
and Az). An asymmetric cut is made on q ■ Ax because 
bremsstrahlung distorts the shower shape in the r~(f> di- 
rection. Since the magnetic field bends an electron's tra- 
jectory, but not a photon's, the bremsstrahlung photons 
tend to enter the calorimeter to the side of the primary 
electron opposite the bending, which is determined by 
the electron's charge. By multiplying the charge by Ax, 
most of the distortion from bremsstrahlung photons is 
isolated to q ■ Ax < 0. 



B. Plug electron identification 

The electron clusters in the plug region, subtending 
1.2 < \ridet\ < 3.0, are limited to 2 x 2 towers (two tow- 
ers in pseudorapidity by two towers in azimuth). Since 
the Moliere radius of a typical electron shower is sig- 
nificantly smaller than the size of the plug EM towers, 
the clusters fully contain electron energies. As with the 
central clusters, plug electron clusters are accepted if 
Ehad/Eem < 12.5%. The major difference between cen- 
tral and plug electrons is the tracking. In the central re- 
gion, the COT tracking is very efficient (99.6%), whereas 
in the plug region the efficiency rapidly falls off as \ridet\ 
increases due to the geometrical acceptance of the COT. 
In this analysis no tracking is used for plug electrons. 
The z position of the collision for the event (zt,ertex) is 
provided by the z position of the central electron's track. 
The electron's shower centroid is determined from a fit 
of the energy distribution among the calorimeter tow- 
ers. The direction of the plug electron is determined by 
extrapolating from Zy^rtex to the shower centroid. The 
unmeasured charge of the plug electron is assumed to be 
the opposite of the central electron. The following vari- 
ables are used to discriminate electrons from hadrons in 
the plug region : (1) the ratio of the hadronic energy 
to the EM energy {Ehad/Eem)', (2) the total transverse 
energy within a radius of 0.4 in AR — Ar^g^j -I- Ac/)^ 
of the cluster centroid, excluding the cluster energy itself 
{E^°)] (3) a comparison of the energy distribution in 
3x3 EM towers around the seed tower to the energy dis- 
tributions from test beam electrons (PE'MxIxs)- The 
selection criteria for plug electrons are summarized in 
Table n 

C. Event selection 

Z/7* — > e^e~ candidate events are required to have 
two electrons with Et > 20 GeV. One of the electrons is 
required to be in the central region and to pass the full set 
of identification cuts (see Table^l. This electron is called 
the central-tight electron. The second loose electron is al- 
lowed to be in either the central or plug region and has 
relaxed identification cuts for higher efficiency (see Ta- 
ble P). Based on these selection criteria, two topologies 
are defined for dielectron events: central-central topol- 
ogy, where one central-tight and one central-loose elec- 
tron are required, and central-plug topology, where one 
central-tight electron and one plug-loose electron are re- 
quired. In the central-central topology, the two electrons 
are required to have opposite charge. In the central-plug 
topology, no charge requirement is made since the plug 
electron's charge is not measured. 

The absolute identification efficiencies are measured 
from the Z/j* e^e~ data where one electron is se- 
lected using the central-tight electron cuts, and the sec- 
ond electron is used as a probe to measure the identi- 
fication efficiencies. The total identification efficiencies 
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Variable 


Central-Tight 


Central-Loose 


Plug-Loose 


Et 


> 20 GeV 


> 20 GeV 


> 20 GeV 


PT 


> 10 GeV 


> 10 GeV 


N/A 


Ehad/ Eem 


< 0.055 
+ 0.00045 * E 


< 0.055 
+ 0.00045 * E 


< 0.05 




< 4 GeV 


< 4 GeV 


< 4 GeV 


E/p 


< 2 

for Et < 50 GeV 


no cut 


N/A 


^ shr 


< 0.2 


no cut 


N/A 




< 10 


no cut 


N/A 


q * Aa:: 


> —3 cm, 
< 1.5 cm 


no cut 


N/A 


\Az\ 


< 3 cm 


no cut 


N/A 




N/A 


N/A 


< 10 



TABLE I: Criteria for electron candidates. Z/"/* — > e e~ 
candidates require at least one central-tight electron and at 
least one additional loose electron (central or plug) in the 
event. 
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FIG. 4: Event electron tdentificatton efficiency, {ef)iD, as a 
function of Alee measured from the Z/'y* e'^e~ simulation. 
The dashed line is for forward events and the solid line is for 
backward events. The dip in efficiency below 90 GeV/(? and 
the differences between forward and backward efficiencies are 
due to radiation effects (see Sec. Mil d) . 



are found to be (83.4 ± 0.8)% for central-tight electrons, 
(94.3 ± 0.5)% for central-loose electrons, and (87 ± 2)% 
for plug-loose electrons. Since we measure the ratios of 
forward and backward events to the total events, only 
the relative difference in efficiencies between forward and 
backward events affects the Afb measurement. The rel- 
ative efficiencies and their dependence on are es- 
timated from the Monte Carlo simulation (see Fig. 0J. 
The difference between the forward and backward event 
efficiencies is largest just below the Z pole. The dip in 
efficiency below the Z pole is due to contamination from 
mis-measured events from the pole. The events are mis- 
measured due to photon radiation which also lowers the 
electron identification efficiency. The forward efficiency 
is lower because the events are more forward at the Z 
pole, resulting in more contamination for forward events. 

Based on these selection criteria, we find 1,892 central- 
central and 3,319 central-plug Z/^* e+e~ candidates. 
The invariant- mass distribution is shown in Fig.Oand the 
number of events in each mass bin is given in Table ^ 



Mass Region 


C-C 


C-P 


40 < Afee < 60 GeV/c^ 


69 


85 


60 < Mee < 70 GeV/c^ 


42 


72 


70 < Alee < 78 GeV/c^ 


48 


119 


78 < Alee < 86 GeV/c^ 


204 


329 


86 < A'lee < 88 GeV/c^ 


151 


299 


88 < Alee < 90 GeV/c^ 


301 


512 


90 < Alee < 92 GeV/c^ 


416 


610 


92 < Alee < 94 GeV/c^ 


330 


543 


94 < Mee < 100 GeV/c^ 


243 


545 


100 < Alee < 105 GeV/c^ 


30 


68 


105 < Alee < 120 GeV/c^ 


29 


61 


120 < Alee < 140 GeV/c^ 


13 


31 


140 < Alee < 200 GeV/c^ 


9 


36 


200 < Alee < 300 GeV/c^ 


6 


8 


300 < Alee < 600 GeV/c^ 


1 


1 



TABLE II: The number of central- central (C-C) and central- 
(C-P) Z/^* — + e^e^ candidates in each mass region. 



D. Detector simulation of Z/^* e^ e 

The simulation is used to understand the detector's ge- 
ometric acceptance for electrons, bremsstrahlung radia- 
tion from interactions with the detector material, and the 
energy scale and resolution for electrons in the tracking 
chambers and calorimeters. In this section, comparisons 
of the kinematic and geometric distributions between the 
data and simulation are discussed, and any discrepancies 
are presented as a possible systematic uncertainty. 

The geometric acceptance is largely dependent on the 
location of the interaction vertex (z^ertea; distribu- 
tion) and the geometric description of the detector. Fig- 



ure |HJ'V shows the distribution of r\det for electrons in the 
Z/7* — > e~ data sample compared with the predic- 
tion from the Monte Carlo simulation. The discrepancies 
shown in Fig. are used to estimate a systematic un- 
certainty on the fiducial acceptance (see Sec. IVIII . The 
Monte Carlo simulation correctly models the observed 
^vertex distribution (Fig. |Sj3). Fitting the Zuertex distri- 
bution to a Gaussian function yields an average position 
in z of -1-2.5 cm and a width of 28 cm. 

The calorimeter energy scale and resolution in the sim- 
ulation are tuned so that the mean and width of the 
Z e^e~ peak in the simulation are consistent with 
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FIG. 5: Distribution of dielectron invariant mass from 
Z/"f* e^e~ candidates in 72 pb^^ of Run II data. 



those from the data (see Sec lIIIEfl . The Et distribution 
of the electrons for central-central and central-plug events 
after the simulation tuning are shown in Figs.EP andEjD. 
(The background prediction included in these plots is dis- 
cussed in Sec. II VI ) 

The amount of material between the collision point 
and the outer cylinder of the COT is tuned so that 
the electron identification variables sensitive to external 
bremsstrahlung in the data match with those in the sim- 
The material in the Monte Carlo 



ulation (see Sec. IIII F(l 
simulation between the interaction point and the track- 
ing volume is tuned using the E/p distribution as shown 
in Fig.[7| The ratio between the number of events in the 
high tail {E/p > 2) and that in the peak is used to cali- 
brate and determine an uncertainty due to the modeling 
of the material in the detector simulation. 



E. Electron energy scale and resolution 

Both local and global energy scale corrections are ap- 
plied to the electron energy. Local corrections are applied 
to improve resolution by correcting for variations in the 
energy response of the calorimeter. They include cor- 
rections for time dependence, variations in the response 
at different points within a calorimeter tower [l^ . and 
variations in the gains of the different calorimeter tower 
channels. Electrons from the W sample and the inclu- 
sive electron sample are used to calibrate these varia- 
tions. The reference for correcting the electron energy 
is the track momentum as measured by the COT. Uni- 
formity is achieved by adjusting the tower energy re- 
sponses (gains) until the mean E/p is flat as a function 
of time and (j>, and agrees with the Monte Carlo simu- 
lation as a function of rj. Figure |H1 shows the invariant- 
mass distributions near the Z peak for central-central and 
central-plug events for data and Monte Carlo simulation. 
High-energy electrons are measured with a resolution of 
£^ = 13|% ® 1.5% [ll| in the central calorimeters and 
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1 1 Z/y^ e*e MC 
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Electron E., (GeV) for C-C events 



Electron E^ (GeV) for C-P events 



FIG. 6: The electron ridet(A) and event vertex z position 
(B) are important for understanding the geometric acceptance 
in the simulation. The electron Et distribution for central- 
central(C), and central-plug(D) are used to the kinematic dis- 
tributions. The Monte Carlo histograms are normalized to 
the number of entries in the background-subtracted data. The 
detector simulation for electrons has been tuned so that the 
Z peak and width match the data (Sec. MU bXi . The discrep- 
ancies are used to estimate systematic uncertainties on the 
Monte Carlo simulation (see Sec. I Wl) . 



• Data 
m Z/y^ e*e- MC 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

E/p for electrons 

FIG. 7: The E/p distribution for central electrons in the 
Z/^* e"^e~ data sample. Points and histograms are data 
and Monte Carlo simulation, respectively. The Monte Carlo 
histogram is normalized to the number of entries in the data. 
The high tail (E/p > 2) of this distribution is used to cali- 
brate the amount of material in the Monte Carlo simulation 
between the interaction point and the tracking volume. 



E 



^ © 0.7% 



[Sfl in the plug calorimeters. 



F. Electron charge identification 

The charge measurement of electrons is essential for 
this analysis because the charge determines the sign of 
cos 9* (see Eq. and therefore determines whether 
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FIG. 9: The rate of same-sign events as a function of the di- 
electron invariant mass from the Z/-y* —^ e^e~ Monte Carlo 
sample. The rate of charge misidentification is half of the 
same-sign event rate. The last bin shows the rate for all 
events with Mee > 300 GeV/c^ . A fit to a constant between 
100 GeV/(? < Mee < 300 GeV/c^ has been made with the 
result that Z/'y* events have an average same-sign rate (pO) of 
1.45%. 



FIG. 8: Invariant-mass distributions ofZ/j* — > e"'"e~ candi- 
dates: two electrons in the central region (top), one electron 
in the central region and the other in the plug region (bottom). 
Points and histograms are data and Monte Carlo simulation, 
respectively. Energy scale corrections and extra energy smear- 
ing are applied to the simulation so that the Gaussian widths 
and peaks match. 



the event is forward or backward. In the central-central 
topology, we measure the charge of both electrons and 
require that they have opposite sign. The central-central 
opposite-sign requirement removes events where one of 
the electrons has a misidentified charge from the event 
sample, leaving very little ambiguity on the forward- 
backward measurement due to charge misidentification. 
However, in events with a central-plug topology, the 
charge of the two electrons is determined solely from the 
central electron. So a misidentification of the central elec- 
tron's charge will switch the sign on cos 0* , and result in 
a misassignment of the event as forward or backward. 
For this reason, the charge misidentification rate needs 
to be properly understood, and well modeled by the sim- 
ulation. 

We use the Z/j* — > e+e^ Monte Carlo sample with 
the central-central topology to study the sources of 
charge misidentification and to measure the misidenti- 
fication rate. The misidentification rate is determined 
by counting the number of events where both electrons 
have the same charge. If the rate is small, the proba- 
bility of having same-sign events is approximately twice 
that of misidentifying the charge of a single electron. Fig- 
ureElshows that above the Z pole {M^e > 100 GeV/c^), 
the rate of events with the same sign is approximately 



fiat up to Mee ^ 300 GeV/c^. The drop in the rate 
at the Z pole and the increase in the rate in the next 
lower bin are attributable to hard bremsstrahlung fol- 
lowed by e'^e~ pair production in the material. This 
process is referred to as charge misidentification due to 
"trident" electrons. Figure ITUl shows a "trident" electron 
where a positron radiates a hard bremsstrahlung pho- 
ton in the material, and the photon subsequently con- 
verts into an electron-positron pair in the material. The 
electron from the photon conversion carries the highest 
momentum, and the charge of the primary electron (e"*") 
is assigned to be negative (e^). The Monte Carlo sam- 
ple shows that the charge misidentification rate coming 
from "trident" electrons is (0.7±0.3)%. The other source 
of charge misidentification is the drift chamber tracking 
resolution of S( — ) ~ 0.001 GeV^^. At higher ener- 

Prp 

gies, the tracks become almost straight (— ^ 0), and 
the charge determination has a higher probability of be- 
ing wrong. The last bin, which includes all events with 
Mee > 300 GeV/c^, has a misidentification rate mea- 
sured from the simulation of (1.1 ± 0.2)%. Corrections 
for charge misidentification are included as part of the ac- 
ceptance calculation via simulation. The dominant sys- 
tematic uncertainty comes from the uncertainty in the 
amount of material between the interaction point and the 
tracking volume. A comparison of the same-sign events 
between data and Monte Carlo simulation is used to get 
an estimate on the background. The comparison shows 
agreement at the Z pole where there is little background 

(Fig.cnj. 
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FIG. 10: A schematic diagram of a "trident" electron where a 
positron radiates a hard bremsstrahlung photon in the material 
and the photon converts into an electron-positron pair in the 
material. The electron from the photon conversion carries the 
highest momentum, thus the charge of the primary electron 
(e~^ ) is assigned to be negative (e^). 



G. Charge dependence of electron efficiencies 




-0.5 



0.5 

for W electrons 



A systematic bias in the forward-backward asymme- 
try may occur if the detector response to electrons differs 
from that to positrons. We compare acceptances and effi- 
ciencies for electrons and positrons using the ^ e^v 
Monte Carlo and data samples. In Fig. the num- 
ber of e"*" events and e~ events are plotted as a function 
of q ■ rjdet of the electrons. The difference between the 
<Z • 'Hdet < region and the q ■ rjdet > region comes from 
the intrinsic charge asymmetry in production in the 
proton direction and the anti-proton direction. Because 
the average momentum of m (u) quarks is larger than that 
of d (J) quarks in the proton (anti-proton), VF+ {W~) 
event production is enhanced in the proton (anti-proton) 
direction. 

Small differences are observed in the data between the 
6+ and e~ yields fFig. ITTB'). These could be caused by 
differences in the detector response between electrons 
and positrons, differences between the \r]det\ > and 
IVdetl < detectors, or an asymmetric Zyertex distribution 
about 0. Similar differences are seen in the Monte Carlo 
sample (Fig. where the differences come purely 

from the asymmetric z^ertex distribution (see Sec. Ill C|l . 
The asymmetry between the e"*" and e~ yields due to 
the asymmetric Zy^rtex distribution is measured to be 
(0.7 ± 0.2)% from the simulation. After the effect of the 
Zvertex offset is taken into account, the observed asymme- 
try in the data is reduced to (0.1 ± 0.7)%. The impact of 
this asymmetry on the forward-backward charge asym- 
metry measurement is therefore estimated to be small 
compared to other sources of systematic uncertainties. 



IV. BACKGROUNDS 

The dominant sources of background to the process 
pp — !■ Z/"f* e+e^ + X are: 

1. Dijet events where the jets are misidentified as elec- 
trons. 



FIG. 11: The q-r]evt distribution of the electrons in W Monte 
Carlo (A) and data (B) events. The asymmetry about is 
due to the charge asymmetry in the W production. The dif- 
ference between and e~ is largely due to the asymmetric 
distribution of the event vertex (zvertex) distribution. 



2. W -\- X ^ ev -\- X , where AT is a photon or a jet 
misidentified as an electron, 

3. Z/^* — > T+T^ — > e+e^i^T-i^epT-i/e, 

4. W+W- -> e+e-UeVe, 

5. W^Z where Z e+e", 

6. ti — > VeVehh. 

The determination of ApB requires knowledge of the 
number of background events and the forward-backward 
charge asymmetry of the background events in each mass 
bin. 



A. Dijet background 

Dijet events are the dominant source of background in 
our sample. This background category includes events 
with two high-py jets and those with one jet and one 
photon. Hadrons in jets can be misidentified as elec- 
trons, and jets can contain real electrons produced in 
semi-leptonic heavy- flavor decays (from c or 6 quarks). 
Photons from initial- or final-state bremsstrahlung are 
identified as electrons if they are in the plug region. In 
the central region, photons must convert in order to be 
identified as electrons because of the tracking require- 
ment. The physics objects producing electron candidates 
in dijet events are discussed in Sec. IIVAII The forward- 
backward charge asymmetry from the dijet background 
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is discussed in Sec. lIV A2l In Sec. lIV A3l the rate of jets 
faking electrons is measured, the dijet background in the 
data is estimated, and the invariant-mass distribution of 
the dijet background is extracted. In Sec. IIV A 41 the 
amount of dijet background in the central-central topol- 
ogy is cross-checked using same-sign events in the Z j^* 
data and Monte Carlo samples. 



1. Sources of electron fakes from dijet Monte Carlo 

The Monte Carlo sample of ~750,000 dijet events is 
used to understand the sources of high-p^ fake electrons. 
We find 47 central-tight electrons, 179 central-loose elec- 
trons, and 1702 plug-loose electrons in this sample. We 
then look for hadron jets or photons with Et > 15 GeV 
that match the direction of these fake electrons. The 
dominant sources for central-tight electrons are photons 
which have converted and light quark jets (those origi- 
nating from u, d, and s quarks). Light quark jets are 
the dominant source for the loose electrons, especially 
the plug electrons. The fraction of events found in each 
category is shown in Table ITTll 



2. Charge correlation and Afb distribution of dijet events 

We do not expect any correlation between the charges 
of the two fake electrons in the dijet events. This expec- 
tation is checked using the hadron-enriched data sam- 
ple, where both electrons are "jet-like"; events collected 
with the high-py central electron triggers that have two 
electromagnetic clusters passing the kinematic cuts for 
electrons in this analysis and a significant amount of en- 
ergy near them {E'!f° > 4 GeV). An additional cut of 
< 10 GeV eliminates a large fraction of the re- 
maining W contamination in the sample. The se- 
lection criteria are summarized in Table IIVI The sam- 
ple contains 8595 (8797) events where the two electron 
candidates have the same (opposite) charge. Although 
the difference is statistically significant {2.2a), there are 
only 2% more opposite-sign events than same-sign events. 
This demonstrates that the charges of the two electrons 
identified in dijet events are nearly uncorrelated. 

The raw forward-backward asymmetry of the dijet 
background is measured with the same hadron-enriched 
sample in the 15 mass regions. As shown in Fig. 1121 
the A^pg values are close to zero, consistent with the 
symmetric angular distribution expected for dijet events. 
This analysis assumes A'^p^ = for the dijet background, 
where the events are split evenly between forward and 
backward categories when subtracted from the candidate 
sample. The measured asymmetries listed in Table IVl are 
used for estimating a systematic uncertainty on this as- 
sumption. 




FIG. 12: 
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A^g of dijet background estimated from the hadron- 
sample. 



3. Estimation and A/ee distribution of the dijet background 

The dijet sample used to calculate the fraction of jets 
faking electrons (the single-jet fake rate) must pass the 
20 GeV single-jet triggers, have two jets with Et > 
20 GeV, have < 10 GeV, and have no more than 
one loose electron in the event. Jets are clustered us- 
ing a cone size of AR = 0.4. These requirements ensure 
that the electroweak contamination from W and Z de- 
cays to electrons is negligible. The fake rate is defined as 
the fraction of jets in the sample which pass the electron- 
selection criteria. The measured rate is plotted as a func- 
tion of jet Et in Fig. El Due to the bigger cluster size 
of jets compared to electrons, jet energies are larger than 
the corresponding fake-electron energies. For example, 
when a jet of Et = 25 GeV fakes an electron, the elec- 
tron candidate is typically reconstructed with an Et of 
20 GeV. The single-jet fake rate with i;^"^* > 25 GeV 
(or E^ > 20 GeV) is measured to be (2.7 ± 0.2) x 10"* 
for central-tight electrons. Note that jets used in this 
measurement are a mixture of both triggered and non- 
triggered jets. 

The single-jet fake rate for non-triggered jets can be 
measured using jets with Et below trigger thresholds. 
We use jets with Et < 45 GeV in the jet 50 GeV data 
sample (50 GeV single jet triggers) and Et < 95 GeV 
in the jet 100 GeV data sample (100 GeV single jet trig- 
gers). In addition, Monte Carlo dijet events without trig- 
ger requirements are studied for comparison. Table IVII 
summarizes the rates in the four different samples. The 
fake rate for jets without trigger biases are roughly one 
half of the fake rate measured in the 20 GeV dijet sam- 
ple. The difference in rates is expected since the trigger 
requirements select jets that deposit a significant frac- 
tion of their energy into one tower and hence look more 
"electron-like." 

The number of dijet background events in the 
Z/7* — !■ e+e" sample is estimated by applying the av- 
erage single-jet fake rate for a tight electron and a loose 
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Dijet Events 



Electron Type 


Light- Quark Jets 


Heavy-Quark Jets 


Gluon Jets 


Photons 


Central- Tight 


0.38 ± 0.09 


0.15 ± 0.06 


0.13 ± 0.05 


0.34 ± 0.09 


Central-Loose 


0.63 ± 0.06 


0.08 ± 0.02 


0.20 ± 0.03 


0.09 ± 0.02 


Plug-Loose 


0.82 ± 0.02 


0.01 ± 0.00 


0.07 ± 0.01 


0.10 ± 0.01 



W ei/ Events 



Electron Type 


Light- Quark Jets 


Heavy-Quark Jets 


Gluon Jets 


Photons 


Central- Tight 


0.28 ± 0.09 


0.00 ± 0.00 


0.00 ± 0.00 


0.72 ± 0.15 


Central-Loose 


0.63 ± 0.09 


0.01 ± 0.01 


0.01 ± 0.01 


0.35 ± 0.07 


Plug-Loose 


0.20 ± 0.03 


0.01 ± 0.01 


0.01 ± 0.01 


0.78 ± 0.07 



TABLE HI: The physics objects producing electron candidates in Monte Carlo dijet and W ev events and their fractional 
contributions for central-tight electrons, central-loose electrons, and plug-loose electrons. The fractional contribution of the 
different physics objects producing electron candidates varies depending on the type of electron and background type. Light- 
quark jets are those originating from u, d, and s quarks, and heavy-quark jets are those originating from to b and c quarks. 
Uncertainties include only Monte Carlo statistics. 



Hadron-enriched selection 



Pj, > 10 GeV/c (central only) 



^corrected ^ GcV 



> 4 GeV 



< 10 GeV 



TABLE IV: Selection criteria for the hadron-enriched data 
sample. 



electron to every pair of jets in the 20 GeV dijet sample 
that pass the kinematic cuts. The systematic uncertainty 
on this estimation is determined by taking the simplified 
case where each event has two and only two jets, one 
triggered and one non-triggered. Let / denote the fake 
rate without trigger biases and 2/ the average fake rate 
in the 20 GeV dijet sample (as described in the previous 
paragraph). Then, the fake rate for the triggered jets is 
3/. The event fake rate is estimated to be 4/^ if the 
average fake rate is applied to both jets (used for the 
estimation of the dijet background), and it is 3/^ if the 
rates for the non-triggered and triggered jets are applied 
separately. The difference between these two estimates, 
which is 25%, is taken as a systematic uncertainty on the 
single-jet fake rate due to the jet trigger bias. 

The invariant-mass distribution of the dijet back- 
ground is determined from the 20 GeV dijet sample 
where at least one of the two jets is found in the cen- 
tral region. Jet energies are corrected to represent elec- 
tron energies where the correction factors are determined 



from Gaussian fits to the 



distributions for each of 



the three categories of electrons. (The Gaussian fits 
to -~t have means ranging from 0.88-0.92, and have 

widths ranging from 0.05-0.09 depending on the type of 
electron.) Figure 1711 shows the invariant-mass distribu- 
tion of the central-central and central-plug dijet events, 



Mass Region 
GeV/c^ 


# Events 


^FB™ (dijet) 


C-C 


C-P 


40 < Mjj < 60 


5.4 ± 1.9 


51 ± 18 


0.00 ± 0.01 


60 < Mjj < 70 


1.5 ±0.5 


25.7 ±9.0 


0.00 ± 0.01 


70 < Mjj < 78 


0.7 ±0.2 


15.3 ±5.4 


-0.01 ± 0.01 


78 < Mjj < 86 


0.4 ±0.1 


11.1 ±3.9 


0.01 ± 0.01 


86 < Mjj < 88 


0.1 ±0.04 


2.2 ±0.8 


-0.09 ± 0.03 


88 < Mjj < 90 


0.1 ±0.04 


1.9 ±0.7 


0.00 ± 0.03 


90 < Mjj < 92 


0.1 ±0.04 


1.8 ±0.6 


-0.02 ± 0.03 


92 < Mjj < 94 


0.1 ±0.04 


1.6 ±0.6 


-0.02 ± 0.03 


94 < Mjj < 100 


0.1 ±0.04 


4.1 ± 1.4 


-0.03 ± 0.02 


100 < Mjj < 105 


0.1 ±0.04 


2.6 ±0.9 


-0.02 ± 0.02 


105 < Mjj < 120 


0.2 ±0.07 


4.8 ± 1.7 


-0.04 ± 0.01 


120 < Mjj < 140 


0.1 ±0.04 


2.9 ± 1.0 


-0.02 ± 0.02 


140 < Mjj < 200 


0.1 ±0.04 


1.9 ±0.7 


-0.05 ± 0.02 


200 < Mjj < 300 


0.0 ±0.00 


0.2 ±0.07 


-0.01 ± 0.03 


300 < Mjj < 600 


0.0 ±0.00 


0.0 ±0.00 


-0.07 ± 0.11 


Total 


9.0 ±3.2 


128 ± 45 


N/A 



TABLE V: The estimated number of events and mea- 
sured A^pg of dijet background in each invariant-mass bin. 
The number of events is estimated using the rate of jets 
faking electrons. The A^pg (dijet) values are measured 
from the hadron-enriched data sample. When calculating 
AFB{Z/-y* — > e'^e~), A^g (dijet) — is assumed. The mea- 
sured A^g (dijet) is used to estimate a systematic uncertainty 
on the measurement. 



weighted by the probability of faking a dielectron candi- 
date event. Table shows the number of dijet events 
expected in each invariant-mass bin. The total num- 
ber of dijet background candidates is estimated to be 
9.0 ± 3.2 for the central-central topology, and 128 ± 45 
for the central-plug topology. The invariant-mass dis- 
tribution of the dijet background along with the Monte 
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Central-Tight 


Central-Loose 


Plug-Loose 




Rate (xlO^*) 


Rate (xlO"") 


Rate (xlO"'') 


20 GeV dijet 


2.7 ±0.2 


13.7 ±0.6 


51.6 ± 1.0 


50 GeV dijet 


1.3 ±0.2 


6.5 ±0.4 


30.2 ± 1.2 


100 GeV dijet 


1.7 ±0.3 


5.5 ±0.5 


27.6 ± 1.8 


Monte Carlo 


1.3 ±0.3 


7.6 ±0.9 


24.8 ±3.1 



TABLE VL The rate at which a jet fakes an electron. The 
jets in the 20 Ge V dijet sample are a mixture of triggered and 
non-triggered jets. The rate in the 50 GeV and 100 GeV dijet 
samples is measured only for jets with Et below the trigger 
threshold, Et < 45 (95) GeV for the 50 (100) GeV dijet 
sample. The rate in Monte Carlo dijet events is measured for 
jets with Et >20 GeV. 



0.0014 E 
0.0012 E 

0.001 E 

0.0008 E 
0.0006 E 
0.0004 E 
0.0002 E 



X^/ndf 17.49/17 
pO 0.0002711 ±2.347e-05 
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+ 



Central Jet E ^ (GeV) 



0.0035 




/ ndf 


29.35 / 20 


0.003 




pO 


0.001366 ± 5.931 e-Q5 
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0.01 
0.008 
0.006 
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22.37 / 22 


pO 


0.005157 ±9.834e-05 



(C) 



Plug Jet E^ (GeV) 



FIG. 13: The rate of jets faking electrons for central-tight (A), 
central-loose (B), and plug-loose (G) requirements in the 20 
GeV dtjet sample. The fake rate is measured for jets with 
Ei^^ > 24 GeV for central-tight and plug-loose electrons, and 
Ei^^ > 26 GeV for central-loose electrons to take into account 
relative energy differences between jet clustering and electron 
clustering. pO is the value returned by the fit for a zeroth order 
polynomial. 



Carlo Z/^* — > e+e prediction is shown in Fig. 1211 



4. Gheck of dijet background using the charge of the 
dielectrons 



We cross-check the estimation of the dijet background 
for the central-central topology using same-sign events 
in the data sample after correcting for the contribution 



Mean 62.27 
RMS 22.16 
Integral 8.978 




Dijet M.. (GeV/c ) for central-central 



> 

(5 



«10 



UJ10 



10 



10 



Mean 69.85 
RMS 23.64 
Integral 127.7 




50 



100 





1 








1 




250 


300 


35C 


400 



Dijet M.. (GeV/c ) for central-plug 



FIG. 14; The invariant-mass distribution measured in the 20 
Ge V dijet sample. The single-jet fake rate is applied to each 
jet m the dijet pair. The jet energy has been corrected to 
represent electron energies (see Sec. IIV A 3l l. 



from charge misidentification. The dominant contribu- 
tion to the charge misidentification of electrons comes 
from "trident" electrons except at very high energies 
where the curvature resolution is the dominant source 
The number of same-sign — > e^e~ 
■'trident" electrons is estimated from the 
The Monte Carlo samples are 



(see Sec. lTTm 
events due to 
Monte Carlo simulation, 
normalized such that the number of opposite-sign events 
in the Z peak is the same between the data and the 
Monte Carlo simulation. Figure \\b\k shows the result- 
ing invariant-mass distributions in data and Monte Carlo 
simulation for opposite-sign events and same-sign events. 
The number of same-sign events is found to be 36 in the 
data and 23 in the simulation, resulting in the estima- 
tion of 13 dijet background events for the central-central 
topology. Figure [T5l 3 shows the invariant-mass distribu- 
tion for the same-sign data events after the Monte Carlo 
same-sign events are subtracted. A systematic uncer- 
tainty of 4 events is estimated by repeating this pro- 
cedure with the Monte Carlo samples generated with 
various amounts of material. Hence the dijet back- 
ground in the central-central topology is estimated to be 
13±6(stat.)±4(syst.) events. This is in good agreement 
with the estimate obtained from the fake-rate method of 
9.0 ± 3.2 events. 
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FIG. 15: (A) invariant-mass distributions of data (points) and 
Monte Carlo events (histogram) in the central- central topology 
where the Monte Carlo sample is normalized to the data using 
opposite-sign events in the Z peak. (B) invariant-mass distri- 
bution of same-sign data in the central- central topology after 
subtracting the same-sign MC distribution. The Monte Carlo 
subtracted same-sign distribution in (B) is a measurement of 
the dijet background. 



B. Electroweak and top backgrounds 

The electroweak and top background events are es- 
timated using the Monte Carlo simulation. Table IVllI 
shows the theoretical cross section and the number of 
events expected in the sample of 72 pb^^ for each pro- 
cess. The systematic uncertainties on these background 
estimates reflect the 6% uncertainty on the integrated 
luminosity, the 5% uncertainty on the acceptances, and 
the uncertainties on the theoretical production cross sec- 
tions. The simulated events for each process which pass 
the selection requirements are used to determine both the 
invariant-mass distributions and the expected forward- 
backward asymmetries. The invariant-mass distributions 
for Z T~^T~ and W-\-X ev-\-X are shown in Fig. 116! 

The dominant electroweak background comes from the 
W+X eiy+X process where X is a photon or hadronic 
jet and the Z/j* —>■ t^t^ process where r — > evr^e- 
While the Z/^* t^t^ e~ Vri^eVT^e background 

is reliably simulated by LO-based generators such as 
PYTHIA, iheW + X ^ ev + X background estimated 
from PYTHIA requires cross-checks with the next-to- 
leading order (NLO) calculations. PYTHIA, using the 
parton-shower algorithms to generate photons and jets, 



Process 


a ■ BR (pb) 


Events Expected 


C-C 


C-P 


Dijct 


N/A 


9.0 ±3.2 


128 ± 45 


W -\- jet/'y ev -\- jet/'y 


2, 690 ± 100 


1.8 ±0.2 


25.4 ±2.4 


Zh* "* r+T- 


252 ±9 


5.6 ±0.5 


7.2 ±0.7 


W+W- e+e'veu. 


0.15 ±0.01 


1.5 ±0.1 


1.8 ±0.2 


W^Z where Z e+e" 


0.15 ±0.01 


1.4±0.1 


1.7 ±0.2 


tt e^vv -\- bb 


0.08 ±0.01 


1.1 ±0.1 


0.7 ±0.1 



TABLE VII: Summary of expected backgrounds. Cross sec- 
tions for the electroweak and top processes are taken from the 
following reference s: rp.pl/ for W and Z, ^ for W+W' 
and W Z , and \3l\. for ti. Monte Carlo estimates are 
normalized to 72 pb~^ . 




Integral 12.8 



(GeV/c ) for Z-^xx 




(GeV/c") for W+X 



FIG. 16: Invariant-mass distributions for (A) Z/'y* t^t 
and (B) W -\- X ev -\- X (X = q,g or y) backgrounds. 



is expected to give lower cross sections than the Stan- 
dard Model prediction for high-p^ photons and jets. We 
cross-check the PYTHIA calculation with two matrix- 
element calculations, ALPGEN for ly + 1 parton pro- 
duction and WGAMMA for T4^ ± 1 photon production. 
The number of W + X ev + X background events 
estimated from PYTHIA is 27 ± 5. The combination of 
ALPGEN and WGAMMA estimates this number to be 
24 ± 3, in good agreement with the PYTHIA expecta- 
tion within the statistical uncertainties. The difference 
between the PYTHIA and the ALPGEN/WGAMMA es- 
timates is used as a systematic uncertainty on Aps- Ta- 
ble shows the sources of the electron candidates in the 
PYTHIA W + X sample, indicating that the dominant 
source is from photons faking an electron (mostly in the 
plug region). 
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V. ACCEPTANCE AND CORRECTIONS 

In order to correct the raw forward-backward asymme- 
try measured in Z/7 decays and obtain A^g'^ , this analy- 
sis must account for any acceptances and detector effects 
that can change the number of forward and backward 
events found in each M^^ bin. Although there is very lit- 
tle in the detector and analysis that treats forward and 
backward events differently, the angular distribution of 
the events, radiation, and detector resolution can each ef- 
fectively change event acceptances and the reconstructed 
invariant mass so that net differences arise for forward 
and backward events. The dominant contributions to the 
detector acceptances are kinematic acceptance (aj^^), ge- 
ometric acceptance (a^go„j), and electron identification 
efficiency (ef^)- The major effects that contribute to 
event migration between Mge bins are the energy reso- 
lution of the detector (af^g), and internal and external 
bremsstrahlung (a^g^^)- Although the simulation is ulti- 
mately used to determine the corresponding corrections, 
generator-level distributions are used to separate and un- 
derstand the various effects fSecs. IV ij| and IV C|l . 
These studies assume the Standard Model ApB distri- 
bution from PYTHIA, and are included to demonstrate 
how the acceptance and detector effects independently 
change A™^. 

We present two unfolding methods to reconstruct 
A^^'^ . In Sec. IV Dl the simulation is used to find a pa- 
rameterization to transform Af^' into A'j?^. This pa- 
rameterization is used to determine the A^ p^'^ t hat cor- 
responds to our measured A^pg. In Sec. IV El correc- 
tion factors are calculated which transform the number 
of forward and backward events in the candidate sam- 
ple to the number of forward and backward events that 
existed prior to being degraded by detector acceptances 
and Mee bin migration. The bias from the input A^g^ on 
the correction factor calculation is also discussed. Both 
unfolding analyses assume the Standard Model -jjj— to 
be able to calculate the event migration effects. 




cos e* (40 < M„ < 78 GeV/c^) COS e" (40 < M^^ < 78 GeVIc') 




cos e* (78 < M^^ < 105 GeVIc') COS 9* (78 < M^^ < 105 GeV/c^) 




cos e* (1 05 < M„ < 200 Ge V/c^) cos e* (1 05 < M^^ < 200 GeV/c^) 



FIG. 17; For three mass ranges, the cos 9* distributions (A, C, 
E) before (open histograms) and after (filled histograms) the 
geometric and kinematic cuts; and the acceptance, t[^eom''Jfci„> 
as a function of cos 9* (B, D, F). 



tributions are asymmetric, due to the forward-backward 
charge asymmetry. When integrating the positive and 
negative portions of cos0*, the convolution of the Drell- 
Yan spectrum with the acceptance gives a different total 
acceptance for forward and backward events. For exam- 
ple, the detector has a low acceptance for events with very 
high |cos 6* I because of the polar coverage, and in the case 
of high-mass events, this removes a greater percentage of 
forward events than backward events (see Fig.llThi'l. The 
Mee dependence of the fiducial and kinematic acceptance 
is shown in Figs.ll8l'V andll8b. where the forward events 
have a higher acceptance than the backward events below 
the Z pole, and vice versa above the Z pole. 



A. Fiducial and kinematic acceptance: ag^om and 
at. 

Kin 

The kinematic and fiducial event requirements sculpt 
the polar angular distribution of the outgoing electrons 
and positrons, especially at high |cos0*|. The distribu- 
tions in Fig. El show this effect for Monte Carlo events 
prior to the simulation of detector and radiation effects. 
The invariant-mass range is split into three different re- 
gions, the low-mass region (40 < Mee < 78 GeV/c^), 
the Z pole region (78 < Mee < 105 GeV/c^), and the 
high-mass region (Mee > 105 GeV/c^). These are the 
regions where ApB is roughly at the low extremum, mid- 
dle point, and high extremum, respectively. Although 
the acceptance for these requirements is nearly symmet- 
ric for positive and negative cos 6* , the initial cos 6* dis- 



B. Corrections for energy resolution: a^es 

The energy resolution of the calorimeter causes mis- 
measurements of the invariant mass which can place 
events in the wrong invariant-mass bin. If the asym- 
metry is changing with invariant mass, events placed in 
the wrong invariant-mass bin alter the measured asym- 
metry in that bin. The effect of the energy resolution on 
the measured asymmetry is therefore largest near the Z 

pole where ^"^J ^ is largest and the bin sizes are smallest. 

As events from the Z pole are smeared out, some events 
move to higher masses where ApB is more forward, and 
some events move to lower masses where the ApB is more 
backward. This event migration increases the correction 
factor for backward events above the Z pole and increases 
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FIG. 18: The Mee dependence of acceptances; a^som (A), a^^i^ 
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(B), a*gs (C), and a^^^ (D). The kinematic acceptance, a^.-^, 
is calculated after geometric cuts are applied, the acceptance 
due to resolution, a^es; calculated after geometric and kine- 
matic cuts are applied, and the acceptance due to radiation, 
'^r^ad' calculated after geometric and kinematic cuts and res- 
olution smearing are applied. 



the correction factor for forward events befow the Z pole. 
The effect of the energy resolution on the correction fac- 
tor is shown Fig.llHC. This plot is made by smearing the 
energy of the generator-level electron according to the 
measured detector response. 



C. Corrections for bremsstrahlung: a^^^ 

The invariant mass can also be mis-measured due to 
the effects of final state bremsstrahlung. Most of these 
photons are emitted co-linearly with the electrons and 
deposit their energy into the same calorimeter towers as 
the electron. In these cases, the energy of the electron is 
recombined with its radiation products. Photons which 
are not recombined with electrons in the energy measure- 
ment lower the measured invariant mass and can cause a 
candidate event to land in the wrong invariant-mass bin. 
The effect is expected to be the most significant just be- 
low the Z pole. Figure I18D demonstrates the effect of 
internal bremsstrahlung on the acceptance at the gener- 
ator level (no external bremsstrahlung is included in this 
plot), when photons are not recombined with electrons. 
The correction factor below the Z pole is large because 
of events from the Z pole that are being mis-measured 
due to radiation. The Aps at the Z pole is larger, and 
boosts the correction factor for forward events. 



D. Parameterized acceptance and smearing 

In order to transform A^'^" into A^pg, a parameter- 
ized function is needed to take into account the event 
loss and migration of events from one invariant-mass bin 
to another. From the Monte Carlo simulation, we com- 
pute the efficiency efj^ for an event that is forward and 
in Mee bin i to be found in the detector as forward in 
bin j. Likewise, 15 x 15 efficiency matrices are calculated 
for backward events being found as backward (e^^), for- 
ward as backward {^fj^), and backward as forward i^fj^)- 
So the diagonal elements of e^^ and e^^ approximately 
represent the acceptance for events in those bins, and the 
off-diagonal elements approximately represent the event 
bin migration. The elements of e^^ and e^^ are much 
smaller and represent those events that are reconstructed 
with a cos 9* of the wrong sign. These efficiencies have a 
small dependence on the underlying cos 6* distribution, 
and so there is some residual dependence on the input 
A^g'' used to calculate them. This might have a small 
bias on the result if the cos 9* distribution in the data is 
very different from the Standard Model. In Sec. IVII Al 
the cos 9* distribution shows very good agreement be- 
tween data and the Standard Model prediction. The four 
matrices are defined as 



FF 



(sel) 



Nf (gen) 
Nf{sel) 
Nfigen) 



BB 



BF 



Nf{sel) 
Nfigen) 
Nfisel) 
Nfigen) 



(16) 



(17) 



where Nfigen) [Nfigen)] is the number of forward 
[backward] events generated in the i-th Mee bin, and 
iNf isel) [Nfisel)] is the number of forward [backward] 
events selected in the z-th Mee bin. By solving for or 
N^ in Eq. Q , we find the expected number of forward 
and backward events for a given Aps, 

Nf^'iphys) = N^iphys) ■ (0.5 ± 0.5 • iAfy%). (18) 

Nl°*iphys) is the total number of events in bin i, and 
is computed from the Standard Model . Then by 
combining the expected number of events and their cor- 
responding efficiencies, the expected number of events 
found in the detector is written as 

Nfiraw) = f;(6,7 . Nfiphys) + eff ■ Nfiphys)ll9) 

i 

Nfiraw) = Y.^ef^ ' + ^l/" ' Nfiphys)l20) 

i 

For a given A or 15 A^^'^ values, Nf (raw) and Nfiraw) 
can be used together with Eq. (0 to find iA'^pg)j, which 
can be directly compared to the A'^pg measured in the 
detector. 



(^Fb). =5(A,j) = 



Nfiraw) - Nfiraw) 
Nfiraw) + Nfiraw) 



(21) 
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E. Calculation of the correction factor af^^ 

The multiplicative correction factor (a+^)i [(a~or)i], 
used to correct ApB '^^ Eq. ifT^ . is designed to be multi- 
plied by the number of forward [backward] events in the 
i-th Mee bin. It corrects for event losses and for events 
that migrated into bin i from another bin with a differ- 
ent ApB- The full simulation is used to calculate the 
which include the effects of fiducial and kinematic accep- 
tance, energy resolution, bremsstrahlung, and electron 
selection efficiency. The combination of all these effects 
is defined as (a^r)i Eq. (j^ . 



N!"'{sel) 



(22) 



In Fig. I19t \ and Table IVIIII the resulting correction 
factor is shown as a function of the dielectron invariant 
mass. This result can be compared to the studies done at 
the generator level for each individual effect (see Fig. ll8|l . 
The overall acceptance is lower due to the more detailed 
fiducial cuts and electron identification efficiencies in the 
full simulation. It is important to note that the cor- 
rection for event migration depends on the input A^^^ 
assumption used to calculate (0^^^)^. The correction can 
only be calculated by assuming the difference in A^g'^ 
between the bin being corrected and the bins from which 
the events migrated. If the A^^" is significantly differ- 
ent in the two bins and there is a large amount of event 
migration, the result of the correction can be biased by 
this assumption. 

To measure the bias of the input Aps on the final re- 

^^Corrcctcd result 

suit, the correlation matrix ( . A'cfcptancc input ) is measured 
by separately reweighting the input ApB ii^ each bin in 
the Monte Carlo sample used to calculate the correction 
factor. The modified correction factors are then used to 
correct the Afb from a high-statistics pseudo-experiment 
using Eq. H15|l . and the change in the ApB result in each 
Mee bin is measured relative to the change in the input 
ApB (Fig. Iiyb). The lowest mass bin and four highest 
mass bins are fairly independent of the A^p^'^ used to cal- 
culate the acceptance, whereas the bins near the Z pole 
have strong correlations. 

To minimize the bias, one of the results uses the A^g'^ 
from Sec. IVII Dl calculated using Standard Model con- 
straints, as the input to the (0^^^)^ calculation. In this 
approach, the uncertainties on the input A^^g'^ from the 
coupling fits are taken as systematic uncertainties on the 
acceptance measurement. 



VI. SYSTEMATIC UNCERTAINTIES IN Afs 

Systematic uncertainties due to energy scale, energy 
resolution, the amount of passive material in the detec- 
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FIG. 19: The correction factors, a^or (j^)- The input 
is taken from the Afb couplings fit fSec. IVII Dl) . The corre- 
lation matrix between the input A'^g" to the correction factor 
calculation and the corresponding change in the result due to 
Eq. H5i (B). The axes are labeled by the Mee bin number 
where bin 1 is the lowest mass bin (40< Mee <60 GeV/(?), 
bin 15 is the highest mass bin (300< Mee <600 GeV/c? ), and 
bin 7 IS the Z pole (90< Mee <92 GeV/c^). 



tor, and background estimation are considered. For a 
given source of uncertainty, a change is made to the cor- 
responding parameter in the simulation, and the impact 
on the asymmetry is evaluated after that change. The 
difference between the asymmetry with the changed pa- 
rameter and the nominal one is taken as the uncertainty 
from that source. The change in the input parameter 
is either a one standard deviation {la) uncertainty on 
the variable in question or a change in an assumption on 
that input. The systematic uncertainties depend on the 
method that is used to extract A^pg^. There are sep- 
arate tables for the systematic uncertainty on A^^^'' for 
the two ApB measurements fTables llXl and lX|l . The Aps 
measurement based on the af^^ correction factors has an 
additional uncertainty due to the Standard Model as- 
sumptions used the calculate these factors (see Table lx)l . 
The systematic uncertainties are quoted with a sign that 
represents the sign of the change in ApB due to the -l-lfj 
variation and is not used in the measurement of ApB- 
A positive variation represents an increase in the energy 
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Mass Range 


Forward 


Backward 


40 < < 60 GeV/c^ 


0.08 ± 0.00 


0.07 ± 0.00 


60 < < 70 GeV/c^ 


0.34 ± 0.05 


0.25 ± 0.01 


70 < Mee < 78 GeV/c^ 


0.53 ± 0.05 


0.36 ± 0.01 


78 < Mee < 86 GeV/c^ 


0.72 ± 0.06 


0.53 ± 0.02 


86 < Mee < 88 GeV/c^ 


0.72 ± 0.03 


0.60 ± 0.02 


88 < < 90 GeV/c^ 


0.39 ± 0.01 


0.35 ± 0.00 


90 < Mee < 92 GeV/c^ 


0.16 ± 0.00 


0.16 ± 0.00 


92 < Mee < 94 GeV/c^ 


0.31 ± 0.00 


0.33 ± 0.00 


94 < Mee < 100 GeV/c^ 


0.56 ± 0.01 


0.65 ± 0.02 


100 < Mee < 105 GeV/c^ 


0.39 ± 0.00 


0.43 ± 0.01 


105 < Mee < 120 GeV/c^ 


0.37 ± 0.00 


0.40 ± 0.00 


120 < Mee < 140 GeV/c^ 


0.39 ± 0.00 


0.40 ± 0.01 


140 < Mee < 200 GeV/c^ 


0.42 ± 0.00 


0.43 ± 0.01 


200 < Mee < 300 GeV/c^ 


0.45 ± 0.00 


0.45 ± 0.01 


300 < Mee < 600 GeV/c^ 


0.42 ± 0.00 


0.43 ± 0.01 



TABLE VIII; The correction factors, a^or, for the different 
mass bins using Afb from the Standard Model coupling fits 
of Sec. WIIIA 




> 
a 



-2-10 1 

Loose Electron ri^ 




-1 1 

Loose Electron ri 



det 



scale and resolution, an increase in the amount of mate- 
rial, and an increase in the amount of background. 

The following systematic uncertainties have been in- 
vestigated and have been found to have negligible ef- 
fects on the measurement of ApB'- fiducial acceptance, 
charge misidentification, ^^^(dijet), and trigger efh- 
ciency. These effects are not included in the total sys- 
tematic uncertainty. 



FIG. 20: The peak (A) and width (B) of a Gaussian fit to the 
Mee peak between 86 and 98 Ge V/(? as a function of rjdet of 
loose electrons. The points are the data, and the lines are the 
chosen variations in the energy scale and resolution used to 
estimate the systematic uncertainties. 



B. Systematic uncertainty from energy resolution 



A. Systematic uncertainty from energy scale 

Variations in the calorimeter energy scale can lead to 
events being placed in the wrong invariant-mass bins. For 
example, near the Z pole where the asymmetry is increas- 
ing monotonically with respect to the invariant mass, a 
positive variation in the energy scale will cause a system- 
atic decrease in the asymmetry (or a systematic increase 
for a negative variation in the energy scale). In general, 
a variation in the energy scale will have an effect only in 
the region where the bins are of the same order as the size 
of the variations (or smaller) and where the asymmetry 
is changing. In this analysis, uncertainties due to the en- 
ergy scale are expected only in the region near the Z pole. 
Figure QUI shows the Gaussian peak of the invariant mass 
as a function of the rjdet of the electron. Based on the 
distribution of calculated masses, the central calorimeter 
scale is varied by 0.5% and the plug calorimeter scale is 
varied by 1% to estimate the systematic uncertainties. 
The chosen energy scale variations are shown as lines in 
Fig. 1201 The corresponding shifts in ^^^^ are shown in 
Tables HXI and 1X1 



Variations in the energy resolution impact the forward- 
backward charge asymmetry in much the same way as 
variations in the energy scale. Instead of systematically 
shifting the events upwards or downwards, they tend to 
smear the forward-backward charge asymmetry to an av- 
erage of the bins around the bin in question. For ex- 
ample, a positive variation in the resolution near the Z 
peak will cause a systematic decrease in the measured 
asymmetry above the Z peak and a systematic increase 
in the asymmetry below the Z peak. As in the case of 
the energy scale, only narrow bins in the region where 

the "^^^^ is large will be affected. Figure EH shows the 
Gaussian width of the invariant mass as a function of 
the rjdet of the electron. Based on the distribution of the 
widths, the central and plug calorimeter resolutions var- 
ied by 0.5 GeV in the central calorimeter and 1.5 GeV in 
the plug calorimeter. In the central calorimeter, only the 
change due to increasing the resolution is used in calcu- 
lating the systematic. The chosen variations are shown 
as lines in Fig. [201 The corresponding shifts Aps are 
shown in Tables IIXI and H 
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C. Systematic uncertainty from amount of material 

The amount of material in front of the calorimeter af- 
fects the energy measurement of electrons. The uncer- 
tainty in the amount of material is estimated to be less 
than 1.5% of a radiation length (Xq) in the region be- 
tween the interaction point and the tracking volume, and 
less than ^Xq in the region between the collision point 
and the plug calorimeter. 

The systematic uncertainty on Afb due to the mate- 
rial is estimated by changing the amount of material in 
the simulation. The changes include adding or subtract- 
ing an extra 1.5%Xq of copper in a cylinder at 34 cm 
(just before the COT) and ■^Xq of steel on the face of 
the plug calorimeter. The corresponding shifts in ApB 
are shown in Tables Hxl andlXl The bins most sensitive to 
the amount of material are those just below the Z pole; 
this is more easily seen in the systematics on the fit result 
(Tabic EI). 



Mass Range 
(GeV/c2) 


Energy 
Scale 


Energy 
Resol. 


Material 


Bgrnd 


Tot 


Stat 


40-60 


-0.003 


0.010 


-0.035 


-0.034 


0.050 


0.127 


60-70 


0.008 


0.008 


0.033 


-0.065 


0.074 


0.185 


70-78 


-0.018 


0.019 


0.040 


-0.017 


0.051 


0.186 


78-86 


0.037 


0.064 


-0.041 


0.001 


0.085 


0.164 


86-88 


-0.080 


0.112 


0.135 


-0.001 


0.193 


0.211 


88-90 


-0.072 


-0.032 


0.043 


-0.000 


0.090 


0.108 


90-92 


0.058 


0.042 


-0.024 


0.000 


0.076 


0.064 


92 - 94 


-0.118 


-0.098 


0.049 


0.000 


0.161 


0.161 


94 - 100 


0.031 


-0.077 


0.013 


0.000 


0.084 


0.168 


100 - 105 


-0.046 


-0.146 


-0.103 


-0.000 


0.185 


0.258 


105 - 120 


-0.001 


0.014 


0.045 


0.014 


0.049 


0.138 


120 - 140 


-0.003 


-0.004 


-0.003 


-0.001 


0.006 


0.185 


140 - 200 


0.002 


0.003 


0.012 


0.004 


0.013 


0.165 


200 - 300 


-0.005 


-0.007 


0.006 


0.026 


0.029 


+0.20 
-0.27 


300 - 600 


-0.004 


0.002 


-0.023 


0.000 


0.024 


+0.00 
-0.51 



D. Systematic uncertainty from background 
subtraction 

The central values of Afb are calculated after sub- 
tracting background events in the forward and backward 
regions separately. The number of background events 
in each Afee bin are estimated in Sec. IIVI The system- 
atic uncertainties due to the background estimates are 
taken as shifts in Afb when the estimated numbers of 
background events are varied by their uncertainties. The 
corresponding shifts Afb are shown in Tables Hxl and Ixl 

VII. RESULTS 

A. Comparisons of uncorrected distributions with 
the Standard Model 

The best way to make a direct comparison of the data 
with the Standard Model is to use the simulation to com- 
pare the uncorrected data with simulated events. The 
background events are included using the predicted dis- 
tributions as described in Sec. IIVI The distributions 
from the signal Monte Carlo simulation are normalized 
to the number of events in the data after subtracting 
the predicted background contribution. Only statistical 
uncertainties have been included for the calculation of 
the comparisons. The uncorrected invariant-mass dis- 
tribution from the data is compared to the signal and 
background predictions in Fig. 1211 Since the energy 
scale and resolution in the simulation has been tuned 
to the data at the Z peak fSec. lIIIE)l . the comparison of 
the Mee lineshape gives a slightly better than expected 
X^/DOF=26.7/35. The comparison is also made for the 
cos 9* distribution (Fig. 122(1 in three mass regions where 
Afb is at extremes; 40 < Mee < 75 GeV/c^ where 
Afb is large and negative (giving a x^/DOF=9.9/9), 



TABLE IX: Sources of systematic uncertainty, the total 
systematic uncertainty (Tot), and the statistical uncertainty 
(Stat) on A'^g" when fitting using the parameterization de- 
scribed in Sec. I VZJI fsee Sec. I VijTI) . When a systematic shift 
is tested in two directions (+la and —la of the variable in 
question), the larger shift is chosen. The sign represents the 
sign of the change in Afb due to the +la variation and is 
not used in the measurement of Afb- 



75 < Mee < 105 GeV/c^ where Afb is small (giv- 
ing xVDOF=44.8/39), and Mee > 105 GeV/c^ where 
Afb is large and positive (giving x^/DOF=7.8/9). Fi- 
nally, the comparison of in 15 Mee bins (Ta- 
ble IXI|) with the Standard Model simulation gives a 
X^/DOF=15.7/15. The data shows excellent agreement 
with the Standard Model in all of these distributions. 
The objective of the following sections will be to obtain 
the corrected A^g'^ and Z couplings that can be used 
without the CDF simulation. 



B. The Standard Model prediction 

Currently, there are a number of programs that gen- 
erate Drell-Yan events produced in hadron collisions. 
PYTHIA generates events using leading-order (LO) cross 
sections and incorporates initial-state QCD radiation and 
initial- and final-state QED radiation via parton shower 
algorithms. HERWIG uses LO cross sections with ini- 
tial state QCD radiation via parton shower algorithms. 
ZGRAD |37| includes full 0{a) electroweak corrections 

but no QCD corrections, resulting in j)^^* ~ 0. The 
gluon resummation program VBP 38, 39], which does 

the gluon resummation in the qt space at low p^^''* 

and reduces to NLO QCD at high p^^^* , does not in- 
clude any electroweak corrections. Unfortunately, there 
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Mass Range 
(GeV/c2) 


Energy 
Scale 


Energy 
Resol. 


Material 


Bgrnd 


Input 


Tot 


Stat 


40-60 


-0.013 


0.013 


-0.028 


-0.053 


-0.018 


0.065 


0.108 


60-70 


0.012 


0.012 


0.027 


-0.038 


0.074 


0.089 


0.095 


70-78 


-0.010 


0.014 


-0.022 


-0.009 


-0.041 


0.050 


0.072 


78-86 


-0.015 


0.033 


0.012 


-0.001 


-0.053 


0.066 


0.043 


86-88 


-0.013 


0.015 


0.013 


0.001 


-0.034 


0.041 


0.048 


88-90 


-0.011 


0.005 


0.011 


-0.000 


-0.014 


0.021 


0.035 


90-92 


0.006 


0.006 


-0.009 


0.001 


0.009 


0.015 


0.031 


92 - 94 


-0.009 


-0.005 


-0.010 


0.000 


-0.004 


0.014 


0.033 


94 - 100 


-0.014 


-0.022 


-0.010 


0.000 


-0.017 


0.033 


0.034 


100 - 105 


-0.017 


-0.080 


-0.045 


0.002 


-0.015 


0.095 


0.099 


105 - 120 


-0.013 


-0.014 


0.027 


0.011 


-0.004 


0.035 


0.091 


120 - 140 


-0.004 


-0.004 


0.008 


0.023 


-0.008 


0.026 


+0. 14 
-0.15 


140 - 200 


0.004 


0.004 


0.011 


0.006 


-0.007 


0.016 


+0.14 
-0.15 


200 - 300 


-0.008 


-0.008 


0.017 


0.033 


-0.005 


0.039 


+0.18 
-0.24 


300 - 600 


-0.017 


0.017 


-0.032 


0.000 


0.000 


0.040 


+0.00 
-0.64 



TABLE X: Summary of uncertainties on A'^g" calculated using Eq. 115\) and the couplings fits for Afb as input to the 
acceptance calculation. When a systematic shift is tested in two directions f+lcr and —la of the variable in question), the larger 
shift is chosen. The sign represents the sign of the change in Afb due to +la variation and is not used in the measurement of 
Afb- The systematic uncertainty due to the Standard Model assumptions used in the correction factor calculation is labeled as 
"Input. " 
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FIG. 21: Invariant-mass distribution of the data compared 
to the prediction for signal and background combined. The 
points are the data, the histogram is the signal Monte Carlo 
sample, and the shaded histograms are the background predic- 
tions. The contributions are added or stacked. 



CTEQ5L parton distribution functions, and PYTHIA 
with MRST2001 parton distribution functions. In each 
mass bin, a band is constructed to extend from the lowest 
to the highest values of the six ApB calculations. The 
uncertainty on the theoretical prediction due to different 
event generators is taken to be the width of the band 
in any mass bin. In Fig. 1231 each calculation is com- 
pared to the center and width of the band in each Mee 
bin. A comparison of these Monte Carlo programs also 
demonstrates the extent to which the CoUin-Soper frame 
(Sec. Ill D|) minimizes the impact of the transverse mo- 
mentum of the incoming quarks. The PYTHIA (LO) and 
ZGRAD programs, which have no initial-state QCD radi- 
ation, can be compared to the VBP, ZGRAD-hPYTHIA, 
and PYTHIA programs which include initial-state QCD 
radiation. The difference in ApB between having and not 
having initial-state QCD radiation is negligible compared 
to the ApB measurement uncertainties (see Fig. I24II . 



is no one program that includes both 0{a) electroweak 
and NLO QCD corrections. A calculation that includes 
0{a) electroweak and some QCD corrections can be ob- 
tained by running ZGRAD with the parton showering 
code in PYTHIA. Six Monte Carlo programs are used 
to constrain the possible values for the Afb measure- 
ment. They are PYTHIA, VBP, ZGRAD, ZGRAD + 
PYTHIA, and PYTHIA with no QCD corrections with 



C. Afb measurement without Standard Model 
constraints 

The forward-backward charge asymmetry is measured 
using the maximum likelihood method and comparing 
A^PB returned by the parameterized acceptance and 
smearing from Sec. IV Dl to the data at detector level. 
This method is unbiased since it does not make any prior 
assumption about the values of ApB. The number of 
events in various invariant-mass bins between 40 GeV/c^ 
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FIG. 22: Distributions of cos 6* for the three mass regions 
of the data compared to the predictions for the signal and 
background combined. The points are the data, the open his- 
tograms are the predictions from signal Monte Carlo simula- 
tion, and the shaded histograms are the predictions from back- 
ground. The contributions are added or stacked. 



and 600 GeV/c^ and A^pg are summarized in Table IXll 
The probability for finding a forward or backward 
event for a given Aps is 



(23) 



and the probability to find forward and N-^ back- 
ward events among N events is given by the binomial 
distribution: 

P{Afb,N^,N^)^ i%.)-{l + AFBf" 

•(l-AFBr^-(i)^^+^^(24) 

Defining a as the negative logarithm of the likelihood 
function, we obtain: 

15 

a = E -Nf ■log{l+g{A,i)) 

i 

■log{l~g{A,i)) + C, (25) 
where A = {Af^B^h, • ■ • , (<r)i5}- (26) 
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FIG. 23: Each of six theoretical calculations compared to the 
center and width of the theoretical band. The variation from 
different theoretical calculations is expressed as the width of 
the band which is determined by the highest and lowest values 
of ApB, including the uncertatnttes, in each mass bin. 



The large correlations between invariant-mass bins 
near the Z pole (see Fig. Et^) lead to large uncertain- 
ties OJi ApB- In this case, it is custoniarv to add a reg- 
ularization function to the likelihood l4(|| to reduce the 
variances. We choose the Tikhonov regularization, which 
adds a function proportional to the second derivative of 
the true distribution. 



ftrueiMee, A^) 



dAL. 



(27) 



'bin i \ ^^^^ee 

where A, = {(<b'')._i, {A^^^, (<r).+i}, (28) 

and where the integral is performed over the entire range 
of bin i in Meg. Since the binning is not uniform, a 
parabolic form is assumed for ftrue{Mee, 



ftrue{Mee,A,) =a + b- M^e + C ■ 



(29) 



The parameters a, b, and c are solved for by integrating 
/ over bins i — 1, i, and i -\- I, and setting the integrals 
equal to A^. S[Ai] can then be simplified to 



S[A,] = 4 • {c{A,yf ■ AMI 



(30) 



This function is only applied to bins 2 through 9, since 
those are the bins with large migration. Adding this to 
the likelihood, a new equation is made 
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15 



^' = E{-^^ • '05(1 + 5(A, z)) - ■ log{l - g(A, z))} + A • ^ ^[A,] + C. 

1=1 i=2 

I 



(31) 
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FIG. 24: Experimental results for Afb with statistical 
and systematic uncertainties (crosses), and theoretical pre- 
dictions based on six independent calculations as described 
in Sec. I VII B\ (bands). The experimental results for Afb 
are measured by fitting to A^pg wtth Tikhonov regularization 
(Sec. I Vll 61) . The bottom figure shows the Z pole region in 
more detail. 



The parameter A is called the regularization parameter, 
and is arbitrary. A very small A will have no effect on 
the fit, and a very large A will dominate the likelihood. 
In this analysis A is chosen such that a' — a + ^. The 
result is not sensitive to changes as large as a factor of 
2 or i in A. The fitting package Minuit is used to 
as a function of A. A^pg'^ obtained after 



minimize a 



Tikhonov regularization is shown in Fig. 1241 and in Ta- 
bic I^H] for different invariant-mass bins. The statistical 
and systematic uncertainties of the measurement and the 
prediction from PYTHIA are also given in Table IXlTl 



D. Study of the Z-quark couplings using the ApB 
measurement 

As seen in Sec. the vector and axial-vector nature 
of the interaction pp — )■ Z/"f* e+e^ renders the Apg 
measurement a direct probe of the relative strengths of 
the vector and axial- vector couplings between the quarks 
and the Z boson. In this section, the effect of varying 
the Z-quark couplings on the forward-backward charge 
asymmetry of electron-positron pairs is investigated. 

The A'^pg measurement at detector level is compared 
to the following theoretical parameterization: 



A''/^{uv,UA, 



dv,dA) 



A raw 
^FB LO 



A raw 
'^FB LO 
_ A raw 
^FB 0{a) 



{uv,ua, dv,dA) 
(SM) 



{SM), 



(32) 



where A™^lo is the predicted ApB at the leading or- 
der and ^™^0(a) is the ApB calculated at the 0{a) 
electroweak corrections level with the ZGRAD genera- 
tor described in Sec. IVll Bl The quark couplings are 
changed at tree level assuming Standard Model couplings 
for the leptons. Both A^pg lq and 0(0,) are obtained 
after smearing the corresponding v4^'^^with the param- 
eterization presented in Sec. IV Dl The parameterized 
A^p^{uv, UAidv ,dA), function of the Z-quark couplings, 
is compared to the measured A^pg , and the best match is 
used to extract the coupling constants. Table IXlTll shows 
the coupling constants obtained after minimization, 
along with the corresponding statistical and systematic 
uncertainties. The divided by the number of degrees 
of freedom (xVDOF) equals 10.40/11. 

The contribution from the different sources of sys- 
tematic uncertainties on the quark coupling constants is 



given in Table IXI VI and the correlation matrix is shown 
in Table IXVI For a given source of uncertainty, the sys- 
tematic error on the coupling is calculated by repeating 
the fit using the measured shifted by its system- 

atic uncertainty. The difference between the fitted cou- 
pling value obtained using the shifted A^pg and the nom- 
inal one is taken as the systematic uncertainty from that 
source. For the PDF uncertainty, we measured the cou- 
plings with four different parton distribution functions: 
MRST99, MRST2001, CTEQ6L and CTEQ5L (default). 
The uncertainty due to the PDF models is determined by 
the largest difference in the fitted coupling values with 
respect to the couplings obtained with CTEQ5L. The 
sign of the systematic uncertainty is defined in the same 
manner as in Table 1x1 Figure 1531 displays the residuals 
compared to the SM ApB from the data (open markers) 
and from the fit (solid line) as a function of the invari- 
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Mass Range 
(GeV/c^) 


Observed Events 


Background 


A raw 
^FB 

zbstat ± sys. 


COS 6* > 


cos 6* <0 


cose* > 


cos 61* < 


40 < Mee < 60 


76 


78 


37.8 ± 9.2 


32.4 ± 9.2 


-0.09 ± 0.11 ± 0.05 


60 < Mee < 70 


46 


68 


19.0 ± 4.6 


17.4 ± 4.6 


-0.31 ± 0.11 ± 0.04 


70 < Mee < 78 


69 


98 


12.4 ± 2.8 


9.9 ± 2.7 


-0.22 ± 0.08 ± 0.01 


78 < Mee < 86 


267 


266 


8.2 ± 2.0 


7.4 ± 2.0 


0.00 ± 0.04 ± 0.00 


86 < Mee < 88 


246 


204 


1.4 ± 0.4 


1.6 ± 0.4 


0.09 ± 0.05 ± 0.00 


88 < Mee < 90 


420 


393 


1.5 ± 0.4 


1.3 ± 0.3 


0.03 ± 0.04 ± 0.00 


90 < Mee < 92 


550 


476 


1.8 ± 0.4 


1.4 ± 0.4 


0.07 ± 0.03 ± 0.00 


92 < Mee < 94 


481 


392 


1.4 ± 0.3 


1.3 ± 0.3 


0.10 ± 0.03 ± 0.00 


94 < Mee < 100 


463 


325 


4.1 ± 0.9 


3.1 ± 0.8 


0.18 ± 0.04 ± 0.00 


100 < Mee < 105 


59 


39 


2.0 ± 0.5 


1.6 ± 0.5 


0.21 ± 0.10 ± 0.00 


105 < Afee < 120 


67 


23 


4.1 ± 0.9 


3.2 ± 0.9 


0.52 ± 0.09 ± 0.01 


120 < Mee < 140 


29 


15 


2.2 ± 0.6 


2.2 ± 0.6 


0.35 ± 0.15 ± 0.02 


140 < Mee < 200 


29 


16 


3.2 ± 0.6 


1.9 ± 0.4 


0.29 ± 0.15 ± 0.01 


200 < Mee < 300 


11 


3 


0.6 ± 0.2 


0.5 ± 0.2 


0.61 ± 0.22 ± 0.03 


300 < Mee < 600 


2 





0.2 ± 0.1 


0.0 ± 0.0 


1.00^.632 ± 0.00 



TABLE XI: Summary of observed events, estimated backgrounds, and A^pg m the dielectron sample from 72 pb ^ of Run II 
data. The systematic uncertainty on A^pg includes only the background subtraction. 



Mass Range 
(GeV/c^) 


< Mee > 

(GeV/c") 


Measured Apb 


PYTHIA ApB 


40 < Mee < 60 


48.2 


-0.11 ± 0.13 ± 0.05 


-0.214 ± 0.003 


60 < Mee < 70 


64.9 


-0.5llo;l7 ± 0.07 


-0.420 ± 0.005 


70 < Mee < 78 


74.3 


-0.45 ± 0.19 ± 0.05 


-0.410 ± 0.005 


78 < Mee < 86 


83.0 


-0.11 ± 0.17 ± 0.09 


-0.214 ± 0.003 


86 < Mee < 88 


87.1 


0.07 ± 0.23 ± 0.19 


-0.079 ± 0.004 


88 < Mee < 90 


89.2 


0.03 ± 0.13 ± 0.09 


-0.001 ± 0.002 


90 < Mee < 92 


91.0 


0.047 ± 0.077 ± 0.076 


0.054 ± 0.001 


92 < Mee < 94 


92.8 


0.15 ± 0.19 ± 0.16 


0.112 ± 0.002 


94 < A/ee < 100 


96.0 


0.35 ± 0.20 ± 0.08 


0.198 ± 0.003 


100 < Mee < 105 


102.2 


-0.02 ± 0.30 ± 0.19 


0.338 ± 0.006 


105 < Mee < 120 


110.7 


0.67 ± 0.15 ± 0.05 


0.454 ± 0.006 


120 < Mee < 140 


128.2 


0.32l^;J^ ± 0.01 


0.554 ± 0.002 


140 < Mee < 200 


161.2 


0.29lo;l7 ± 0.01 


0.598 ± 0.002 


200 < Mee < 300 


233.6 


0.65lo;27 ± 0.03 


0.609 ± 0.004 


300 < Mee < 600 


352.4 


1 nn+O'Oo+o 

l.UU_Q -0.02 


0.616 ± 0.007 



TABLE XIL Experimental results for Afb measured by fitting to A^pg with a Tikhonov regularization function. Statistical and 
systematic uncertainties are included along with predictions from PYTHIA using CTEQ5L. The uncertainties in the PYTHIA 
predictions are MC statistical errors. < Alee > is the cross section weighted average of the invariant mass in each bin. 



ant mass. The dashed curve corresponds to the residuals 
with the uv quark couphng shifted by la. 

The CDF sensitivity, with ~72 pb~^ of analyzed data, 
is limited, but the values of the couplings are consistent 
with the Standard Model. Figure 1^ shows the contours 
at 68% and 90% confidence level for the u (left) and d 
(right) quark couplings to the Z boson in the vector-, 



axial-vector basis. The closed markers correspond to the 
Standard Model predictions. 

Table IX VII summarizes the effective left- and right- 
handed coupling constants obtained from this study and 
compares them to the Standard Model prediction as 
well as the current experimental values of the effective 
couplings determined from 'model independent' fits to 
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Coupling 


Uncertainty 


Statistical 


Systematic 


Uv 


0.399 


+0.152 
-0.188 


0.066 


UA 


0.441 


+0.207 
-0.173 


0.067 


dv 


-0.226 


+0.635 
-0.290 


0.090 


dA 


-0.016 


+0.346 
-0.536 


0.091 



TABLE XIII; Vector and axial-vector quark couplings with 
statistical and systematic errors obtained from the fit to the 
Afb measurement at detector level. 
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FIG. 25: Residuals compared to the SM Afb as a function of 
the invariant mass from the data ( open marker) and from the 
fit (solid line). The dashed line represents the residuals with 
the Uv coupling from the fit shifted by la. 



neutral-current data [43|. As the present study is dom- 
inated by the statistical uncertainties, the sensitivity of 
these measurements will improve with higher integrated 
luminosities. 

The Z-quark couplings returned by the fit are subse- 
quently used to compute the A^''^^ which is shown for 
different invariant-mass bins in Table IXVIII The uncer- 





En. scale 


Resol. 


material 


bckgrd. 


PDF 


Uv 


-0.056 


-0.023 


-0.025 


0.001 


-0.001 


UA 


-0.013 


-0.048 


-0.009 


-0.044 


0.003 


dv 


-0.013 


0.038 


-0.076 


-0.021 


0.017 


dA 


-0.059 


0.062 


0.006 


-0.025 


-0.018 



TABLE XIV: Contribution of the different sources of system- 
atic uncertainty on the quark coupling constants. 





nv 


UA 


dv 


dA 




1.000 


0.454 


0.303 


-0.037 


UA 


0.454 


1.000 


0.214 


0.428 


dv 


0.303 


0.214 


1.000 


0.548 


dA 


-0.037 


0.428 


0.548 


1.000 



TABLE XV: The correlation matrix of the statistical errors 
in the measurement of the up and down quark couplings. 





This study 


Exp. values [42] 


SM prediction [42] 


Ul 


0.419 


0.330 ± 0.016 


0.3459 ± 0.0002 


dL 


-0.116 


-0.439 ± 0.011 


-0.4291 ± 0.0002 


Ufl 


0.020 tl:{tl 


-0 176 +"°^^ 


-0.1550 ± 0.0001 


dR 


0.105 


-0 023 +'' °™ 


0.0776 



TABLE XVI: Effective left- and right-handed coupling con- 
stants obtained from this study and compared to current exper- 
imental values and the Standard Model prediction from U&]. 



tainties in the resulting A^'^^ are the total (statistical 
and systematic) uncertainties from the fit of the Z cou- 
plings. A^pg^ is then used to calculate the acceptance as 
discussed in Sec. IV El 



Bin 


Mass Range 


Fitted Af^" 





40 < Mee < 60 GeV/c^ 


-0.170 ± 0.074 


1 


60 < Me, < 70 GeV/c^ 


-0.355 ± 0.125 


2 


70 < Mee < 78 GeV/c^ 


-0.373 ± 0.109 


3 


78 < Mec < 86 GeV/c^ 


-0.183 ± 0.124 


4 


86 < Mee < 88 GeV/c^ 


-0.044 ± 0.083 


5 


88 < Mee < 90 GeV/c^ 


0.028 ± 0.053 


6 


90 < Mee < 92 GeV/c^ 


0.088 ± 0.027 


7 


92 < Mee < 94 GeV/c^ 


0.140 ± 0.039 


8 


94 < Mee < 100 GeV/c^ 


0.223 ± 0.060 


9 


100 < Mee < 105 GeV/c^ 


0.342 ± 0.087 


10 


105 < Mee < 120 GeV/c^ 


0.429 ± 0.101 


11 


120 < Mee < 140 GeV/c^ 


0.493 ± 0.135 


12 


140 < Mee < 200 GeV/c^ 


0.506 ± 0.159 


13 


200 < Mee < 300 GeV/c^ 


0.501 ± 0.171 


14 


300 < Mee < 600 GeV/c^ 


0.498 ± 0.175 



TABLE XVII: Forward-Backward asymmetry calculated with 
the Z-quark coupling values returned by the fit. The uncer- 
tainties in Afb are based on the total (statistical and system- 
atic) uncertainties on the Z coupling values returned from the 
fit. 



Assuming Standard Model Z-quark couplings, we can 
also determine the vector and axial-vector couplings be- 
tween the electron and the .Z-boson, using the same 
method. Table IXVf f f I shows the fitted electron coupling 
values along with statistical and systematic uncertain- 
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FIG. 27: Experimental results for Afb along with statistical 
and systematic uncertainties (crosses), and theoretical predic- 
tions based on the six calculations as described in Sec. IVllBI 
(bands). The measured Afb values are corrected for accep- 
tance, efficiency, resolution, and bremsstrahlung. The correc- 
tion assumes the Standard Model and uses the derived A'^g" 



from the Z-quark coupling values, 
the Z pole region in more detail. 



The bottom figure shows 



FIG. 26: 68% and 90% confidence level contours for the u 
(top) and d (bottom) quark couplings to the Z boson. The 
Standard Model predictions for these couplings are indicated 
by closed markers. 



with a xVDOF equals to 12.50/14. The present CDF 
sensitivity on sin^ 9'^^^ is provided by the Z-electron cou- 
plings and is expected to improve with higher statistics. 



ties. Combined LEP and SLD data '4^ as well as the 
Standard Model prediction _42] are shown for compari- 
son. The xVDOF of the fit equals to 13.14/13. For the 
present measurement, the uncertainties are dominated 
by the statistical uncertainties. The contributions of the 
different systematic uncertainties are given in Table lXiXl 
Note that the correlation coefficient between the vector 
and axial- vector couplings is 0.78. 

Finally a fit where the quark and electron couplings 
to the Z boson are expressed as a function of sin^ 9'^^ 
gives: 

sin^ 9%f = 0.2238 ± 0.0040(stat) ± 0.0030(syst), (33) 



E. Afb measurement assuming the fitted Standard 
Model couplings in the acceptance calculation 

The total acceptance, a^^' calculated for each bin 
using the A^g'' obtained from the Z coupling fits (Ta- 
ble The ApB measurements are corrected for ac- 
ceptance, efficiency, resolution and bremsstrahlung using 
Eq. 115|) . The measured Aps and the prediction from 
PYTHIA using CTEQ5L parton distribution functions 
are listed in Table IXXl and the measurements are com- 
pared with the Standard Model theoretical calculations 
(see Sec. IVII B|) in Fig. ^\ This technique, which is bi- 
ased by the Standard Model input, is similar to the Run 
I analysis 
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coupling 


Stat, 
err. 


Syst. 
err. 


Total 
err. 


SLD-f LEP meas. 


SM prediction 


ev 


-0.058 


0.016 


0.007 


0.017 


-0.03816 ±0.00047 


-0.0397 ±0.0003 




-0.528 


0.123 


0.059 


0.136 


-0.50111 ±0.00035 


-0.5064 ±0.0001 



TABLE XVIII: Vector and axial-vector electron couphngs with statistical and systematics uncertainties obtained from a fit to 
the Afb measurement. Combined LEP and SLD data 'Aj] as well as the standard model prediction are also given. 





En. scale 


Resol. 


material 


bckgrd. 


PDF 




-0.005 


0.005 


0.001 


-0.001 


0.001 


CA 


-0.002 


0.055 


0.006 


-0.005 


0.019 



TABLE XIX: Contribution of the different sources of system- 
atic uncertainty on the electron coupling constants. 



physics. The second principal result is a measurement 
of three sets of parameters: the Z-quark couplings, the 
Z-electron couplings, and siv^Ow. All three couplings 
measurements yield results consistent with the Standard 
Model. It may be possible to improve our understanding 
of the Z-quark couplings with a much larger dataset. 
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